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LOAD PLOW CONTINGENCY EVALUATION AND OPTIMAL POWER DISPATCH 
IN LARGE SCALE POWER S YSTEMS -DEVELOPMENT OP NEW AND MORE 

EPPICIENT METHODS 

This thesis concerns itself with the development of new 
and more efficient solution methods for load flow, contingency 
evaluation and optimal power dispatch in large scale power 
systems. Load flow is essential for power system planning, 
operation and control and provides voltages and power flows 
at different points in power system network for different 
operating conditions. Contingency evaluation involves 
determining the security of a system with reference to a set 
of next contingencies and constitutes the main function of 
the power system security analysis. Security functions like 
security monitoring and security analysis are now incorporated 
into the computer program to deal with different operating 
conditions as well as disturbances. On-line monitoring and 
real-time control of the system using process computer is 
essential in order to prevent or minimize costly system 



breakdowns. This led to the development of powerful cal- 
culation procedures for load flow and contingency evalua- 
tion. These techniques must he simple, efficient and fast 
for implementation in real-time on moderate size process 
control computers. The Newton-Raphson method (NM) is the 
basic practical method best suited for off-line studies 
because of its quadratic convergence rate and better accura- 
cy. Numerous other methods of load flow and contingency 
evaluation have been tested and published for off-line and/or 
on-line studies. The fast decoupled technique of Stott in ? 
polar coordinates (FILP) gained popularity over other methods 
due to its fast computing speed and moderate memory require- 
ments for off and on-line calculations. Recently, second 
order load flow techniques have been presented to improve 
accuracy and to reduce computer time of NM. 

The problem of optimal load flow (OLF) is very important 
because of the high fuel and other operating costs. The solution 
of OIF gives the values of control variables such as real 
and reactive power generations and transformer tap ratios for 
which the system operating cost is minimum. The OIF calcu- 
lation is a complex nonlinear programming problem due to the 
large number of variables and constraints. To solve OIF 
efficiently for large practical systems, fast and accurate 
methods of load flow solution are needed. A comprehensive 
comparative analysis of different optimization techniques for 



the solution of OLE was given, by Happ and Sasson. Dommel 
and Tinney proposed a practical and efficient method for 
the calculation of OLE by using KM. Later on, application 
of EDLE in the place of KM is also reported. 

The function of the economic dispatch program is to 
allocate the real power generation among the available on-line 
generating units so that the operating cost is minimized. To 
implement on real-time basis, it requires advances in the 
accurate calculation of penalty factors or a power flow 
solution which is both fast and accurate. This problem has 
been solved with approximate loss formulas as well as more 
accurate and faster methods. 

However, there is still a scope of improving these well 
known methods of load flow, contingency evaluation, optimal 
load flow and economic dispatch in terms of computation time. 

Present day ijower systems are large because of greater 
interconnections. To analyse such large scale systems with 
limited core storage, application of sparsity oriented or- 
dered elimination techniques and/or the various piecewise, 
decomposition and diakoptic techniques is needed. Recently, 
the hierarchical control of power system has been and conti- 
nues to be studied intensively. In this method, the original 
system is decomposed into several smaller subsystems or 
areas coordinated by a central agency. It requires the use 



of small multi-computers with significant reduction in 
computer storage requirements. Sometimes, a net saving in 
computation time is also achieved. A survey of literature 
reveals that the hierarchical solution methods for Dll' are 
limited to the optimization of real power flow only. (Control 
of reactive power flow and transformer tap ratios is negle- 
cted). 

Keeping the above in view, this thesis presents solu- 
tion techniques to the different problems considered. The 
scope of thesis is broadly divided into two parts. 

i) Development of new load flow solution methods along 
with a theoretical study of the convergence characteristics 
for off-line and/or on-line studies and its adaptation to 
contingency evaluation and optimal load flow. 

ii) Development of load flow and optimal load flow methods 
for large scale power systems by a two level computation 
based on hierarchical control theory. 

The chapterwise summary of the work covered in this 
thesis is given as follows : 

Chapter I introduces the problems of load flow, contingency 
evaluation and optimal load dispatch. Then, the theory of 
hierarchical control is presented. A brief review of the 
state of the art and the scope of the thesis are enumerated 
at the end. 



^hapter^II is devoted to the development of the following 
new load flow solution methods, i) East decoupled load flow 
in rectangular coordinates (EDER): The rectangular version is 
derived by linear transformation from EDEP as well as direct 
from EM. It offers computational advantages of 5 to 10 per- 
cent speed over the polar version, while the memory require- 
ments are similar, ii) Newt on -Richard son method (NRM) : 

This method is a variant of EM, obtained by combining both 
the chord method and the Eewt on-iterative method. It is 
about 2 times faster than RM for large systems and has the 
same solution accuracy and convergence rate as EM. The 
Jacobian is readily available for important applications like 
optimal power dispatch, sensitivity analysis and security 
constrained optimization, iii) Network-reduction method (EETR) : 
The method is based on eliminating P-Q buses using sparse 
matrix techniques and applying RRM to an equivalent 
system consisting of only P-V buses. It is comparable to 
EDEP with regard to computing speed and memory space for the 
medium size systems. Test results of the IEEE 14 bus, 30 bus, 
and 57 bus test systems and a heavily loaded 23 bus system 
using PDLR, REM and IEIR algorithms are presented and com- 
pared with the results of EM and EDEP. In addition, 
the results obtained using EDLR for a 153 bus system 
are compared with the result of EDEP. 



Chapter III is concerned with contingency evaluation needed 
for power system security assessment. EDLR and KETR load 
flow algorithms developed in Chapter II are extended for 
the rapid and efficient simulation of network and/or power 
outages. EDXR incorporates the matrix inversion lemma and EETR 
adapts both the matrix inversion lemma and matrix updating algo 
rithm in order to make use of the base case factors to obtain 
all contingency states as in prevalent competing methods 
due to Stott and Peterson. The proposed EDLR algorithm is 
faster than PDIiP by about 10 percent. But the EETR algorithm 
has computational advantages over PDIP if the number of 
P-V buses is roughly about 10 percent of the total number 
of buses In the system. Digital simulation of a number of 
outage cases is carried out on the IEEE 14 bus, 30 bus and 
57 bus test systems., The significant results of these 
system studies are presented and compared with the result of 
EDLP. 

Chapter_IY deals with optimal load flow solutions and eco- 
nomic load dispatch. A new optimal load flow calculation 
method using NRM load flow algorithm and a new fast economic 
load dispatch solution using EETR load flow algorithm in the 
solution approach of Dommel and Tinney are proposed. The 
method of solution using ERM gives an overall advantage in 
speed of approximately 2 to 2£ times compared to the approach 
of Dommel and Tinney without affecting the solution accuracy. 



convergence properties and generality of formulation. 

^he solution approach hy EETR is based on the assumption 
that the voltage magnitude variations at the P-Q buses are 
not too sensitive to the real-power changes and the reactive 
power equations depend weakly on the phase angles of the 
voltages. This makes the algorithm very fast and operate 
with less storage. The optimal load flow and economic dis- 
patch results of the adapted IEEE 14 bus, 30 bus and 57 bus 
test systems using the proposed NRM and RETR algorithms 
respectively are presented and compared with the result of 

Pommel and TLnney's approach. 

/ 

C ha pter V describes a two level computational procedure 
for load flow and optimal load flow calculations in large 
scale power systems based on hierarchical control theory. 
It can be regarded as an application of decomposition. It 
is suitable for multicomputer configuration with less core 
storage requirement. It combines both the advantages of 
sparsity and decomposition. A two level load flow calculation 
method is described using HM, NRM and EDEP. Then, the multi- 
area optimization problem is formulated as an additively 
separable nonlinear programming problem in such a way that 
the optimization of both real and reactive power and trans- 
former tap ratios is possible. The total lagrangian function 
is decomposed into several sub-Eagrangians by introducing 
pseudo variables. The second level determines the pseudo 



variables and the associated Lagrangian multipliers. The 
first level determines the subsystems solutions and the 
associated Lagrangian multipliers. The solution algorithm 
for optimal load flow calculations considering both equality 
and inequality constraints using the overall 
objective function is given in detail. The 
cability of the proposed two level technique is d 
by obtaining the load flow and optimal load flow 
adapted IEEE 30 bus and 37 bus test systems which 
into 2 or 3 subsystems or areas. 

Chapter, VI gives a theoretical study of the conv 
characteristics of the load flow solution methods (ELLR, NEM 
and EETR) developed in Chapter II. The different load flow 
algorithms are expressed into the general nonlinear iterative 
scheme. Eor these algorithms, the convergence conditions, 
which guarantee the existence and uniqueness of load flow 
solution in a specified region of interest, are obtained by an 
application of the convergence theorem used by Wu. These condi- 
tions are verified by using an adapted IEEE 14 bus test 
system for different R/X ratio of the heavily loaded line. 

Chapter VII reviews the major contributions in the thesis. 
Ttirther possibilities of research in this area are also 




indicated 



CHAPTER I 


INTRODUCTION 


1.1 LOAD PLOW STUDIES 

The solution of the load flow problem is carried out 
extensively on a digital computer for power system planning, 
operation and control. This is essential both for off-line 
and on-line applications. The objective of this study is to 
determine the phase angle and reactive power on each P-V bus 
and the phase angle and voltage magnitude on each P-Q bus 
subject to the constraints on the real and reactive power at 
the P-Q buses and the real power and voltage magnitude at 
the P-V buses. The digital solution of load flows has been 
studied extensively using different forms 13,5,8,14,20,901 
of nonlinear equations, The formulation in the nodal admi- 
ttance reference frame is preferred, since it preserves 
sparsity of the network. An excellent review by Stott {25I 
gives the salient features and the comparative merits of 
different load flow solution methods. The important proper- 
ties required of a load flow solution method are high computa- 
tional speed, low storagf, reliable convergence and versa- 
tality. Further, it must be suitable for modern applications 
like security analysis, optimal load flow and economic dispatch 
calculations. 
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1.2 POWER SYSTEM SECURITY 

Recently, a new strategy of security control has emerged 
for the reliable operation of electric power systems and to 
prevent black-outs and blown-outs. This has made a signifi- 
cant impact on power industries and the security functions are 
being implemented with the aid of on-line computers in more 
than 100 power system control centres [ 41 I throughout the 
world . 

1.2.1 Security Concept 

DyLiacco [261 identified three modes for power system 
operation namely the normal state, the emergency state and 
the restorative state. A power system is subjected to two 
sets of constraints, i) load constraints (load flow equa- 
tions) ii) operating constraints (maximum or minimum limits 
on system variables such as voltages, real and reactive 
power generations, transformer tap ratios, phase angle 
differences etc.). By definition, normal state is one in 

t 

which both the load and operating constraints are satisfied. 

A system is in emergency state when the operating constraints 
are not completely met. A system is in restorative state, 
if the operating constraints are satisfied but not the load 
constraints. The normal state may be either secure or in- 
secure. A normal state is said to be secure, if it undergoes 
any contingency without getting into an emergency condition. 

If there is a possible contingency which causes the system to 
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enter the emergency state, then the normal system is in- 
secure, Fig. (1.1) gives a typical scheme for computer 
software organization [41] that links the various security 
oriented functions together. The objective of security 
functions is to keep the system operating in normal state. 

1.2.2 Security Analysis 

Security analysis comprises of two functions : 

(1) Contingency evaluation: It determines the security of 
the system with reference to a set of next contingencies, i.e. 
verifies whether the normal system is secure or insecure, 

(2) Preventive control : It determines suitable preventive 
action, when the system is insecure. 

1.3 0PTIMA1 LOAD FLO¥ STUDIES 

The optimal load flow study determines an operating 
state or load flow solution in which a certain objective 
function like fuel cost or system loss is minimized taking 
the generator real powers, voltage magnitudes and transfor- 
mer tap ratios as control variables subject to the operating 
constraints of the power system. A number of problem for- 
mulations [ 39,70,73 ] exist using different choices of 
objective functions and constraints. In its more general 
form, the problem is a complex non-linear and constrained 
mathematical programming problem. A variety of simplifica- 
tions by linear or quadratic programming have been reported 
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in detail [73,78,106,109 1 . Classical economic dispatch 
can "be considered as a particular case of optimal load flow 
with only generator real powers as control variables and no 
security constraints. Many power system control centres 
today still operate with the above concept, the calculations 
being made on-line every few minutes. 

1,4 HIERARCHICAL CONTROL THEORY 

Hierarchical control methods are becoming popular for 
large scale power system problems. There are two kinds of 
hierarchy which are of great importance. 

(1) Multilayer Hierarchy : This concerns the decomposition 
of the function of system organization into layers. Each 
layer supervises the system with a different model. This 
kind of hierarchy is used in off-line as well as real-time 
applications using on-line centralized computer. 

(2) Multiechelon or Multilevel Hierarchy : In this, the 
large system is decomposed into mpiy subsystems on the basis 
of the physical structure of the system. The overall solution 
is obtained by proper coordination. In general, different 
methodologies presented by Mesarovic [89 land Haimes [104] 

are highly beneficial in the coordination of the subsystems 
of a given decomposition. This kind of hierarchy is useful 
when a decentralized multicomputer configuration (Fig. 1,2) 
is employed with the help of regional or area computers. 




R6'1.2 MULTIAREA TWO LEVEL COMPUTE 
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coordinated "by a central dispatch computer. The decentra- 
lized hierarchical structure has been considered by several 

authors for state estimation [99,107], reactive power and 

I 

voltage control [98,105], economic dispatch calculation 
[87,96,108] and transient stability [100]. 

This thesis addresses itself to the important problems 
of power system just discussed namely, load flow, contingency 
evaluation, optimal load flow and economic dispatch. It 
begins with the . development of new load flow solution methods 
which are competitive and superior in some respects to the 
existing methods. Then, these methods are adapted to con- 
tingency evaluation and optimal power dispatch. Before de- 
lineating the theoretical study of the convergence characteri- 
stics of the new load flow solution methods at the end, a two 
level computational technique based on hierarchical approach 
is developed for load flow and optimal load flow calculations 
of large scale power systems. In the development of solu- 
tion methods for optimal load flow and economic dispatch, 
nonlinear programming techniques are used. 

A brief literature survey of the various models and 
techniques pertinent to these areas is now presented. 

1.5 STATE OE THE ABT 

1,5.1 load Plow 


load flow solution methods are well documented in the 
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literature [ 23,25,40 ] and can "be classified into four cate- 
gories, 

1 ) Iterative methods employing Y^g - matrix or - matrix 

2) Newton methods, 3) Decoupled methods, 4) Other methods. 

The Gauss Seidel method [9 ] using Y BBg - matrix offers 
advantages of simple problem formulation retaining sparsity. 
But it requires a large number of iterations for convergence 
and hence is slow computationally. Further, it presents 
convergence difficulties if the bus admittance matrix is 
not diagonally dominant, which may happen when the line 
configurations have capacitance. The method employing Z-g^g - 
matrix [1,6] with Gauss iterative scheme converges more re- 
liably compared to Y-g Ug - matrix method. Since Z BU g - matrix 
is full, sparsity is lost. Further, additional time is 
needed to generate Z-g US - matrix. To overcome this diffi- 
culty, various diakoptic, decomposition and piecewise methods 
primarily using Z BUg [4,6,14,27] have been developed. 


The Newt on- Raph son method (NM) [2] together with spar- 
sity and ordered elimination techniques [3] is a powerful 
method for obtaining load flow solutions. It uses the first 
two terms of the Taylor series expansion. The convergence 
of this method is fast and reliable. However, the main dis- 
advantage of this method is its large computational require- 
ments due to repeated calculation and triangulation of the 
Jacobian. Meisel and Barnard [12] gave the convergence 
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characteristics of the NM method and proposed the use of 
constant Jacobian for the load flow solution. The algorithm 
requires a larger number of iterations compared to NM and it 
may fail to converge for large systems. The methods proposed 
by Wallach et al [ 15 ] use rectangular coordinates and update 
the inverse of the Jacobian at eveiy kth iteration by appli- 
cation of linear algebra which are not amenable to sparse 

' 

matrix techniques. Several oth er variations of M have also 
been reported [ 23,25,391. Of late, there has been some 

. 

interest in the second order load flow techniques. Sachdev 
and Medicherla [31J extended KM in polar coordinates with 
different variations to include upto the first three terms 
of the Taylor series expansion. This second order approach 
gives better accuracy and also there is some saving in com- 
puting time compared to Ml. Iwamota and Tamura [40] used a 
constant Jacobian NM method in rectangular coordinates based 
on full Taylor series expansion. In this method of formula- 
tion, terms upto second order only are included in the Taylor 
series expansion of load flow equations, the higher order 
terms being zero. .Although this method is claimed to be 
much faster than KM, it requires more memory space. Further, 
convergence difficulties have been reported by the authors 
using flat voltage start. 

After 1970, decoupled methods have received attention for 
the fast a.c. load flow solution. These methods which solve 
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the MW- 6 and WA R-V problems separately, require approximately 
40$ less storage than NM and reduced computer time. Despotovic 
et al [16] developed a decoupled method in -which direct solu- 
tions for the voltage magnitude and phase angle at the system 
buses are obtained. This method was later improved [281 with 
certain approximations. The matrix relating to the calcula- 
tion of voltage magnitude correction is constant and the ma- 
trix relating to the calculation of phase angle correction 
needs re-calculation and re-tri angulation at every iteration. 

A simultaneous-displacements iteration scheme is used which 
reduces the convergence characteristics of the method. Stott 
[19 ] presented a load flow method based on the decoupling 
principle, in which the polar current mismatch equation is 
utilised instead of the conventional WAR mismatch equation 
to improve its convergence rate. It requires the calculation 
of coupling matrices at every Iteration. Making further 
assumptions based on certain network properties, Stott and 
Alsac [291 developed a fast decoupled load flow technique. 

The final matrices are constant and hence the table of factors 
are to be obtained only once at the beginning of the iterative 
process. This method is now acknowledged to be superior to MM 
in terms of computational speed and storage. However, its con- 
vergence rate and accuracy are not better than NM« Recently, 
Berg et al [36] presented a method based on decoupling prin- 
ciples and successive approximation technique. .Although its 
solution ao curacy is similar to that of MM, the resulting 
matrices are not constant. 
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A variety of other load flow solution methods have also 
heen reported. Wallach [7] and Sasson [lOlused minimization 
techniques for load flow solution. These methods, which mini- 
mize an unconstrained scalar objective function usually equal 
to the sum of the squares of the current or power mismatches, 
are not best suited for routine load flow calculations. 
Dusonchet et al [171 proposed a method in which the load flow 
equations are divided into two groups. The first group corres- 
ponding to P-Q buses is solved by the generalized point Jacobi 
method, the second group corresponding to P-V buses by NM. 

This method achieves savings in storage and computational re- 
quirements compared to M. Okamura et al [30] extended EM to 
include the effects of system control devices such as governors, 
exciters etc. Besides, flat and freouency bias tie line con- 
trols have also been discussed. Elangovan et al [ 34 ] presen- 
ted an approximate load flow method based on reduction tech- 

f 

nique without exploiting sparsity. The method, however, does ’ 
not take into account the G-limit violations and its accuracy 
is good only if phase angles at the system buses are small. 

1.5.2 Contingency Evaluation 

The existing methods for contingency evaluation use sensi- 
tivity factors or load flows based on D.C. model, linearized 
A.C. model either full or decoupled and the constant matrix 
model. MacArthur [47 ] and Limmer [49] reported the use of 
sensitivity factors. The linear sensitivity factors known as 
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the distribution factors are generally derived from Xg^g - 
matrix. These factors cannot track the changing topological 
structure of the network and require periodic updating depend- 
ing upon the operating state of the system. Moreover, compu- 
tational and accuracy problems arise. D.C. load flow [511 is 
extremely fast and is comparatively better than distribution 
factors in accuracy. However, it gives only the solution of 
real powers and ignores actual voltage profile in assessing 
over loads. Daniel and Ohen [55 ] developed a method similar 
to d.c. load flow using the changes in real power injection 
at the system buses to simulate a line outage. It does not 
provide any information regarding the reactive power flow. 
Contingency evaluation methods using 1.0. load flow "are. availa- 
ble inrZ^g'as well formulations. 


Stagg and El-Abiad [ 46] investigated the use of Z-matrix 
for contingency evaluation which requires .an excessive core / 
storage. Brown [ 50,60 ] used Z-matrix axes discarding tech- 
nique in order to conserve computer memory and to increase 
the speed of matrix formation for the determination of over- 
loaded lines and low voltage caused by contingencies. These 
techniques are fast and suitable for simultaneous multiple 
outages but do not often yield accurate or complete informa- 
tion regarding the operating state of the power system. 

The Gauss-Seidel method [46] can treat the changes in 
conf iguration very easily but its convergence is very poor. 
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The results obtained using the nonlinear load flow equations 
with KM [54] are very accurate but the computer time is 
excessive. Uemura [53] used constant matrices which are 
obtained by inverting the decoupled Jacobian matrix. Since, 
this method is not amenable to sparse matrix techniques, it 
needs excessive core storage and considerable computer time. 
Stagg and Phadke [52] presented an approach which uses the 
constant Jacobian NM method. The accuracy and speed of the 
method acre not adequate, Sachdev and Ibrahim used sensitivity 
analysis in polar [58] as well as rectangular [61] formula- 
tions. These methods are fast and compute all the contingency 
states using base case factors of the Jacobian. However, some 
additional off-line computations for the suitable set of power 
infection modification factors are also needed. Pai and Khan 
[59] introduced a faster method based on sensitivity approach 
and linear models. Since the resulting matrices are non- 
symmetric, the method requires an increased memory storage. 
Tinney [ 56 ] proposed the use of a compensation method with 
triangular factorization to obtain sparse network solutions 
involving network outages. This led to the development of the 
well known competing methods due to Peterson [571 and Stott 
[29]. The linearized power flow model proposed by Peterson 
et al is based on the direct solution of a sparse system of 
linear network equations. Although this method is fast, it 
gives only an approximate solution. A separate A.O. load flow 
routine is required to obtain a full a.c. load flow solution. 
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Stott and Alsac proposed a fast decoupled algorithm using 
sparse triangular factors of the constant matrices. This is 
essentially a fast method and provides approximate as well 
as full a.c. load flow solutions. 

1.5.3 Optimal load Plow 

A comparative study [ 67] and a comprehensive survey [73 ] 
of different optimization techniques for the solution of op- 
timal load flow (01P) have heen made hy Sasson et al and Happ 
respectively. To start with, Carpentier [64] formulated the 
general OLP problem which tried to meet the optimality rela- 
tionships of Kuhn and Tucker. The method is slow and presents 
a convergence problem. Sasson [66] developed a sequential 
unconstrained minimization technique which presents convergence^ 
difficulties for large systems. Pommel and Tinney [651 pro- 
posed the use of lagrangian multipliers and penalties to asso- 
ciate respectively the equality and inequality constraints to 
the objective function. The method uses NM for obtaining the 
load flow solution. The resulting reduced gradient is then 
used to find corrections for control vector. It is basically 
a first order method and is very sensitive to the selection of 
step lengths used at each gradient step. ALsac and Stott [ 72 ] 
used PDIiP in the place of NM without cheugiug the generality of 
formulation of Pommel and Tinney's approach to reduce memory 
storage and computer time. Sasson et al [68] used penalty 
functions for both equality and inequality constraints which are 
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added "to the objective function. The augmented objective 
function is then minimized by using Hessian method posse- 
ssing second order characteristics to update all the voltage 
magnitude and phase angle of buses. The convergence of this 
method is dependent upon the proper choice of values for the 
penalty factors. Further, the Hessian matrix is large and 
less sparse than the Jacobian matrix in M and its structure 
is quite different from - matrix. Therefore, complex 

programming subroutines are needed to determine its confi- 
guration so that sparse matrix techniques can be employed. 
Rashed and Kelley [71J employed a second order approach to 
find corrections for the control vector. This is an improved 
version of Dommel and Tinney’s method. Billinton and Sachdeva 
[691 used a sequential optimization approach in which the vol- 
tage magnitude is first optimized based on an assumed value of 
real power generation and then using the just optimized vol- 
tage variables, the real power generation is optimized. This 
iterative procedure is repeated until the improvements in the 
optimized variables are reduced below a certain tolerance 
level. Recently, Bala and Thanikachalam [74] proposed a 
method retaining the second order characteristics. An equi- 
valent smaller sized matrix obtained from the original large 
sized Hessian matrix is used to update the control variables. 
Although this method is independent of the proper selection of 
step lengths, the computational requirements per iteration are 
more than in Dommel and Tinney's approach. 
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1.5.4 Economic Dispatch 

Economic dispatch calculationsusing loss coefficients 
(B-coef ficients ) are in prevalent use for real time applica- 
tions. Various techniques [75-77, 79] have been reported 
for the calculation of B-coeff icient s. Shoults et al [84] 
recently presented a method for the computation of B- 
coefficients based upon the method of least squares. It is 
accurate and efficient for solving economic dispatch problem 
for large interconnect ed systems. However, the above prac- 
tice requires a re -computation of loss coefficients, when- 
ever there is a significant change in the network. Alvarado 
[83] proposed a method for the calculation of incremental 
losses and penalty factors using the transposed sparsity 
oriented Newton-Raphson load flow Jacobian. This method is 
claimed to be faster compared to classical B-coefficients 
for large systems. Wadhwa and Handa [81] presented a method 
in which economic coordination equations are solved by NM 
load flow technique to achieve better accuracy and speed, 
Alsac and Stott [80] used EDLP to compute penalty factors 
avoiding explicit Jacobian matrix. Malik and Rao [82] pro- 
posed the use of EDEP in the place of NM to solve the coor- 
dination equations. These methods, which employ EDIiP, reduce 
the core storage and computer time, making it attractive for 
real-time applications. 
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1.5.5 Hierarchical Computation 

There are few methods available in the literature per- 
taining to the application of hierarchical computation for 
solving the load flow and optimal load dispatch problems. 

Happ and Undrill [87J illustrated the basic features of multi- 
computer power flow computation using d.c. power flow solu- 
tion. It is reported that this approach can also be extended 
to an a.c. load flow. Mitter and Dyliacco [851 proposed a two 
level optimization procedure which utilizes sequential un- 
constrained minimization technique to solve the first level 
optimization problem. Happ [86,96] explained how to dispatch 
real power economically through multicomputer configuration. 

This method is an extension of loss coefficients formulation 
in a multi-area case. Another paper by Happ f 97 J describes 
the multilevel piecewise theory for large network problems 
in a generalized manner. Recently, Irving et al [ 107 ] have 
presented their computational experience with decomposed state , 
estimation and optimal real power dispatch based on hierarchi- 
cal control theory. In this method, different linear progra- 
mming techniques have been used for optimal real power control. 
It is reported that the solution times are realistic for on- 
line applications. Bielli et al [108] developed a method for 
real power economic dispatch based on the hierarchical approach. 
This method utilizes quadratic programming at the first level 
optimization and the coordination is achieved at the second 
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level 'by gradient type algorithms. In addition, the opti- 
mization of nondiff erent ial functions and stochastic optimi- 
zation have also been discussed, ill these methods, however, 
consider only .real powers for optimization. The control of 
other variables such as reactive powers and transformer tap 
ratios is ignored. Paranjothi [103] proposed a real and 
reactive power optimization procedure amenable to multi- 
computer configuration. It is reported by the author that 
the method presents convergence difficulties. 

1.6 SCOPE AMD OUTLINE OE THE THESIS 

In this section, the main objectives of the thesis and 
a chapterwise summaiy are given : 

The main objectives of the thesis are as follows : 

i) Development of fast decoupled load flow technique in 

p>i 6 ** ucmM* th, 

rectangular coordinates (EDLR), along the lines ofCFDLPJ 

A 

to improve the computing speed and its adaptation to 
contingency evaluation. 

ii) Development of Newt on-Ri chard son load flow solution 
method (NKM) - a powerful alternative method to NM and 
its adaptation to optimal load flow. 

iii) Development of load flow solution method based on network 
reduction and NRM techniques (NETR), competitive to EDIP, 
in terms of computer time and memory storage for medium 
size systems and its adaptation to contingency evaluation 
and fast economic load dispatch. 
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iv ) Development of two-level computational technique for 
large scale power system load flow and optimization 
based on hierarchical control concepts. 

v) Theoretical study of the convergence characteristics of 
the new load flow solution methods. 

A brief summary of the work reported in the thesis is 
given below : 

In Chapter II, new load flow solution methods namely 
FDLR, ITEM and 1TETR are developed. 

FDLR is derived in two ways. Firstly, through linear 
transformation from FD1P and secondly direct from NM load 
flow equations in rectangular coordinates. The proposed 
algorithm is faster than FDEP for moderately accurate solu- 
tions. 

ITEM is obtained by combining the Chord method [121 and 
the Newton iterative method [13 J. This is an attractive al- 
ternative method for systems where NM is preferred. The 
solution accuracy and convergence rate are practically the 
same as in NM, while the computing time is about 50$ for 
large systems. It requires less storage and number of ite- 
rations for convergence compared to the recently developed 
method [40] in rectangular coordinates without truncation 
error. 

NETR is developed using sparse matrix techniques based on 
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eliminating the P-Q buses and then applying the fast decoup- 
led technique. The reduced system consisting of only P-Y 
buses is solved by NRM. The method is fast, accurate and 
has good convergence properties. The solution accuracy is 
practically the same as in EM. A comparison with EDEP 
indicates both the competitive nature of the method and its 
suitability for on-line applications. The effectiveness of 
the proposed load flow algorithms is investigated on the- IEEE 
14 bus, 30 bus, 57 bus test systems and also on a heavily 
loaded 28 bus system. Besides, EDLR is also tested on a 133 
bus system. A comparison of accuracy, speed and core 
storage requirements of these algorithms is made with EDIiP 
and EM. 

Chapter III deals with contingency evaluation using 
NETR and EDLR load flow algorithms developed in the pre- 
vious chapter. The proposed techniques are very well suited 
for fast on or off-line outage security checks. As in 
prevalent competing methods 1 29 , 57 li these techniques avoid 
the re -tri angulation of the base case matrices for the si- 
mulation of different outage cases. 

In NETR algorithm, the network outages are simulated 
using the matrix inversion lemma and the matrix updating 
algorithm. The technique is sufficiently accurate and is 
comparable to or faster than EDLP depending upon the number 
of P-V buses in the system. In EDLR algorithm, the network 
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outages are simulated with the help of the matrix inversion, 
lemma. This rectangular version is faster and superior in 
accuracy for line flows greater than 50 MW compared to polar 
version i.e. PDLP. 

Digital simulation of a number of network outage cases is 
carried out on the IEEE 14 bus, 30 bus and 57 bus test systems. 
The significant results of these system studies are presen- 
ted and discussed in detail. 

Chapter IV is devoted to the application of NRM and NETR 
load flow algorithms developed in Chapter II for optimal load 
flow and economic dispatch calculations respectively using 
the solution approach of Dommel and Tinney. The method of 
solution using ITEM is faster in speed by more than 2 times 
compared to Dommel and Tinney's approach, without affecting 
the generality of formulation and convergence properties. 

The solution approach by NETR is based on two assumptions, 

1) the voltage magnitude variations of the P-Q buses are not 
too sensitive to the real power changes, 2) the reactive power 
equations depend weakly on the phase angles of the voltages. 
This makes the algorithm competitive in terms of speed and 
core storage for real-time applications. Details of these 
solution algorithms are given. The proposed algorithms are 
tested on the adapted IEEE 14 bus, 30 bus and 57 bus test 
systems. A comparison of accuracy, speed and memory require- 
ments of these proposed algorithms is made with Dommel and 
T inn ey ' s app ro ach . 
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Chapter V concerns Itself with the development of load 
flow and optimal load flow methods for large scale power 
systems based on hierarchical concepts. A two level technique 
utilizing the formulation proposed by Haimes [104] is used. 
Additional pseudo variables are introduced to decouple the 
large scale system. The 1st level determines the subsystem 
■solutions. The second level coordinates via pseudo variables. 
The computational advantages of sparse matrix techniques are 
used in the analysis of subsystems. Firstly, the multiarea 
load flow equations are given and solved by using EM, NRM, 
and EDLP. Then, the multiarea optimization problem is for- 
mulated as an additively separable non-linear programming 
problem by considering the real power generations, voltage 
magnitudes and transformer tap ratios as control variables, 
using the solution approach of Dommel and Tinney [651. The 
effect of acceleration factor and the effect of number of 
subsystems over the rate of convergence are illustrated. 

Details of the optimal load flow solution algorithm consi- 
dering both the equality and inequality constraints using 
the overall cost as the objective function are given. The 
effectiveness of the proposed two level technique is demon- 
strated by obtaining load flow and optimal load flow results 
for the adapted IEEE 30 bus and 57 bus test systems which are 
torn into 2 or 3 subsystems. A comparison of these results 
with those of the original untorn system is also made in terms 
of solution accuracy, storage space and overall computer time. 
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Chapter VI presents a theoretical study of the conver- 
gence charact eristics of the new load flow solution techniques 
(EDER, NHM, NETR) proposed in Chapter II. Each one of these 
load flow techniques is expressed suitably into the form of 
a general nonlinear iterative scheme. Then, the convergence 
conditions for these techniques are obtained by an application 
of the convergence theorem, used by Wu [331 . The validity 
of these conditions is checked on the IEEE 14 bus test syptem. 

Chapter VII highlights the major contributions of the 
thesis, in assessment of the scope for further work is then 


made 



CHAPTER II 


LOAD PLOW SOLUTIONS 

2.1 INTRODUCTION 

Many excellent methods [23,25,40] have been tested and 
published during the last 12 years to solve the problem of 
load flows for large scale power systems. A brief review and 
comparative study of the existing methods have been made in 
Section 1.5.1 of the preceding chapter. In this chapter, new 
load flow solution techniques, competitive to the prevalent 
methods are presented. The salient features of the proposed 
techniques are summarized as given below. 

(1) Past r decoupled load-flow in ^rectangular coordinates (PDLR.), 

PDLR is developed via PDLP using linear transformation and 
also direct from NM load flow equations in rectangular coordi- 
nates. The final matrices turn out to be the same as in PDLP. 

A PDLR iteration is faster than that of a PDLP iteration, since 
trignometric functions do not have to be calculated. The 
existing PDLP program can be easily modified into PDLR. A 
comparison with PDLP indicates that PDLR is faster for mode- 
rately accurate solutions. 

( 2 ) L o ad^ _ f low_ _ s olu ti. on _ by_ Newton-Mchardson method_ ( M] 

NRM is a variant of NM based on the Chord and Newton- 
iterative techniques. The method offers an overall' saving in 
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computing time of approximately 40 to 50^ as compared to EM, 
The solution accuracy and the convergence rate are practically 
the same as in EM. The Jacobian is readily available for the 
important applications like optimal load dispatch, sensitivity 
analysis and security constrained optimization. The load flow 
program can be easily obtained from the existing EM program 
by adding a subroutine for the evaluation of minor iterates. 

(5) Lqad_flqw solution by ^network reduct ;ion and Newton- 

Ri chard son techniques ( EETR) 

HETR is based on eliminating the P-Q buses and applying 
ERM to the P-V buses only, exploiting sparsity of the 
network. The method is comparable to EDLP in terms of computer 
time for medium size systems while the solution accuracy and 
convergence rate are similar to EM. It is very well suited 
for fast on-or off-line outage security checks and economic 
dispatch. 

The proposed algorithms have been tested for the standard 
IEEE 14 bus, 30 bus and 57 bus systems, and also for a heavily 
loaded 28 bus system. PEER algorithm has also been tested for 
a 133 bus system. 

2.2 MATHEMATICAL EO EMULATION OE THE PROBLEM 

The basic problem is to determine the complex bus voltages 
in an A.C. transmission system for a prescribed set of con- 
straints. Table 2.1 gives the classification of buses with 

It becomes necessary to choose one of 


pertinent information. 
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the available P-V buses as slack bus 

to take care 

of the 

system losses. 




TABLE 2.1 CLASSIFICATION OF BUSES 


'Type of 

Specified quantities > 

Unspecified quantities ’ 

bus 



4 


Polar version 

; Rectangular [ 


:: 


version 

P-Q bus 

1. real power consumed ; 

1 .voltage 

; } 

rl.real part ' 

: or 

,2. reactive power i 

magnitude 

: of the : 

'load bus 

consumed 

2. volt age 

voltage 



phase angle 

2. imaginary - 




part of the! 

Mi.:.:. . m- . jnu r tu i 

_ ; 


: voltage ; 

P-V bus 

1. voltage magnitude '■ 

1 . reactive 

; 1. reactive » 

or 

■;2. real power 

power 

; power ] 

voltage con- 

->3. limits of reactive - 

2. volt age 

2. real part 3 

trolled bus 

power 

phase 

5 of the 


> 

angle 

: voltage ? 

' 3. imaginary •; 

* 

f 

\ \ 


part of thef 
l voltage 

t 

1 Slack bus 

’For polar version: ? 

1 .real power 

: l.real power- 

or 

; Voltage magnitude and j 

2. reactive 

! 2. reactive 

1 swing bus 

; phase angle 

power 

power 

or 

■For rectangular version: ; 


* 

■ reference 

Real and imaginary parts . 


< » 

-bus 

of the voltage 


‘ i 

2.2.1 Bus constraint equations 



For 

the ith P-Q bus 




pSP _ p^P _ p - o 
^G-i i Li - - u 


(2.1) 
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Q li - «ti - «1 = 0 (2-2) 

where 

n 

p i + 3 Q i - E i l *i - E i 1 Y ti =5 

D=1 

For the ith P-V bus : 


P S T? _ 

ui 


-*i - ° 


(2.3) 

•y^p _ 

( e i 

* fl 2 )* - 

v i 

(2.4) 

In concise 

vector form. 

the equations (2.1 . 

- 2.4) can be re- 

presented i 

as 




f g( 

[x] , 

)) = 

0 

12. 5) 


where [x] = vector of all unspecified quantities 
[ y ] = vector of all specified quantities 

The load flow problem is therefore reduced to the numerical 
solution of the constraint equations (2.5) to find [x] for a 
given [y] . This basic problem gives only an unadjusted solu- 
tion. However, an adjusted solution can be obtained by sol- 
ving (2.5), subject to the reactive power, voltage magnitude 
and tap ratio constraints. The overall problem is non-linear 
in nature. Hence, iterative procedures are used to obtain 
the numerical solution. 
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2.3 DEVELOPMENT OP LOAD PLOW MODELS 
2.3,1 Newt on-Raphson Method (NM) 

In this basic method, the solution of the non-linear load 
flow equation (2.5) is obtained by linearizing it using the 
Taylor series expansion in the neighbourhood of the initial 
estimates, neglecting the second and higher order terms. The 
polar or rectangular coordinates may be used for the bus vol- 
tages. 

Polar Version : 

The set of linearized equations in polar coordinates 
is represented in a matrix form as 



where 


AP? = real power mismatch for the ith P-V or P-Q bus 


41 - - v i X (( hj 003 + B ij sin VV 


(2.7) 
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AQ? a reactive power mismatch for the ith P-Q "bus 


- «& - v i * « e i 3 si “ e ij ' 


B i3 003 W 
( 2 . 8 ) 


A0 = "bus voltage phase angle, magnitude corrections 

respectively. 


[ «J?1 = submatrices of the Jacobian [J^J 

1 i=1 to 4 

n v = number of P-V buses in the system 
n^ = number of P-Q buses in the system 

and the superscript p refers to the polar coordinates. The 
elements of the Jacobian [ J^ ] are defined in the Appendix A. 

The linear equation (2.6) is solved for the corrections 


t A0 ] 

and 

[AV] to be applied to 

[6 ] and l V] . Then, 

16] 

k+1 

= i e ] k + t A6 ] k 

(2.9) 

[V 1 

k+1 

= [V] k + [ AV] k 

(2.10) 


The iterative calculations are repeated until the real and 
reactive power mismatches are within a specified convergence 


tolerance. 



Rectangular Version : 


The set of linearized equations in rectangular coordina- 
tes ig given in a matrix form as 


[g([X] , 

m) r i - 

U r I [ A x r ] 



(2.1 

r— -j 

; A P r : 

(n v +n L ) 

Dffj.-jux 

i 

1 

f 

— I ~7 T- — . 

_ r 3? r ,M . i 

[ 3f X ? 

.! A Q r • „ 

(n L ) 


t 

_ t 

L J 4l _ 

fcf l ■ a e 

x : 

• A 

4 

K) 


" l 

1 

1 

f J 6 J 

;,*? ; * ? 
~ 3 f J , ; { 




* 


1 W -X 



(n v +n I ) 


(n v 

+n L ) 


where 


Ap? = p^? 
l Gi 


p Sp 

r Li 


L (e.e. + f.f.)G + (f e -e,f )B. . (2.12) 

i-J X J 1 J XJ X J X J XJ 


(for ith P-V and P-Q bus) 


n 


4 <? - 


o^P 

^i 


* (f i e 3 - e iV s i 3 - ' 2 - 13 > 


(for ith P-Q bus) 


2 

At\ = square of the bus voltage magnitude corrections 

(for ith P-V bus) 

and the superscript r refers to the rectangular coordinates. 

r 2? i 

The elements of the Jacobian submatrices (J- 1 are 

i=1,2,...6 

given in the Appendix A. 

liquations (2.6) or (2.11) have been solved for the correc- 
tion vector by exploiting the sparsity • .of' network admittance 
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and Jacobian matrices, through triangular factorization, compact 
storage and optimal renumbering of buses [3J which gives fast 
direct solutions at eveiy iteration. However, MM is computa- 
tionally Inefficient, since it requires the factorization of the 
highly sparse [J] at every iteration. A number of modifications 
[ 25 ] including the recent second order models [31, 40l are avai- 
lable in the literature to reduce the volume of computations. 

2.3.2 Decoupled Methods 

In all the decoupled methods [19,28,29,36], the load flow 
equations have been derived from NM to reduce the memory space 
and computing effort. These methods are based on neglecting the 
coupling submatrices [J^] and [J^J in equation (2.6) on the 
assumption that the coupling between the real power versus bus 
voltage magnitude and the reactive power versus bus voltage 
phase angle are relatively weak;. Among the existing decoupled 
load flow methods, FDLP is simple, reliable and fast and is best 
suited for off or on-line applications. Its precise algorithmic 
form has been determined by extensive numerical studies. A theore- 
tical study of its convergence characteristics has been reported 
later by Wu [33]. 

Past decoupled load flow in polar coordinates (PDIiP) 

The final FDLP equations are given by 

[^] = [B'J [A9] (2. 14) 

l ] * [B n ] [A V] 


(2.15) 
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■where (B'l and f B” 1 are matrices of order (n Tr + n T ) and 
n^ respectively and can be obtained directly from the elements 
of t -B 1 with the following modifications C 29 1 

i) omitting from fB'J the representation of network 
elements namely shunt reactances and off-nominal in-phase 
transformer taps. 

ii) omitting from tB”J the angle shifting effects of phase 
shifters. 

iii) neglecting series resistances in calculating the elements 
of [B’l which then becomes the d.c, approximation load 
flow matrix. 

The real and sparse matrices [B'J and [B") are constant. 
Hence, the table of factors are obtained only once at the 
beginning of the iterative process. Using these constant 
factors, very fast repeat solutions for { A9 ] and[AVJ are 
obtained. The iterative procedure is repeated till the real 
and reactive power mismatches become small enough. Various 
storage schemes can be applied to take the advantage of spar- 
sity and symmetry of [B’l and [B 11 ] . However, [B’l is 
symmetrical only if the phase shifters are neglected. 

2.4 BAST DECOUPLED LOAD PLOW IN RECTANGULAR COORDINATES (EDLR) 

2.4.1 Derivation of EDLR from EDLP by linear transformation 

The real and imaginary parts of the ith bus voltage in 
terms of its magnitude and angle are 
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e i = V ± cos e ; f. = V. sin 0 ± 


(2.16) 


By taking the exact differentials of the equation (2.16), we 
get. 


; 4 e i : i - v i «*» e i ■ 


cos 0^ 


? A0 


: Af il } Y i oos 9 i ! sin e i I 


(2.17) 


Inverting (2.17), we obtain 


. A0 . * 
1 


: ay.! 
; 1 


sin 

_ ^ — 
i 


cos 0^ 


cos 0^ 


sin 6 . * 
i ; 


A e i } 


“i : 


( 2 . 18 ) 


Applying the transformation (2.18) in (2.14) and (2.15) and 
changing CAP] and [ A Q] to rectangular coordinates, the 
following equations are obtained 


[A? r ] 

1 ¥ 1 


CB' 1 C A' : Ag ] tjfl 


(2.19) 


c A 9f 1 

1 V 1 


l A? ! 


I ^ ^ 


! Ae 


(2.20) 


where [ A} ] , [ A£ ] , [A!j| and [ A£ ] are diagonal matrices 

with ii-th element given by 
sin 0 . 

(A^)ii = - y ■ 
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(AJ>) 


COS Oj; 

'2 J 1± = Tf" 


(^Pii = G0S 0 i 
and ^ A 2^ii = sin 9 i 

In the polar version, the voltage magnitude equations 
do not appear for the P-V buses. Hence, the equations (2.19) 
and (2.20) have to be augmented by 


A e 


t AV J = fAij” j A"> ] l £yr 


( 2 . 21 ) 


where [Alt’] and t A^' ] are diagonal matrices with ii-th 
element given by 

(A?' ) 1JL - 2 ; <1? ) u - 2 fj. 

On combining the equations (2.19 - 2.21) and making the same 
assumptions as made in PDLP , namely - 1.0 and 0^ - - 0 

and after some manipulations, we get 


AP, 


; A P r * 
n j 


A 

$ 

U Q r J 

;ih] 

> 

[1] > 

*** ~ 4 

• Ae, 
< 

' • 

j * i = 

• 

it mi 

i N] 5 

' Ae 

2 f 

A Vo ,, 

t S 3 

* 

! [T 1 ' 

I ? 

r 

; Af, 

i • 

i 

* 

A V 2 
m 



»<H 

< 




( 2 . 22 ) 
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^'"'^2.2^2^ *•* ! 


2m nr * 


tH] =: 


*** ^ 


nm nr 


’(B^ 2 e 2 ) ... (B^e m ) | 


^ B 2,m+1 f m+1^ * * * ^ B 2n f n^ 


(~ B n,m+1 f m+1 ^ 
(B 2,m+1 e m+1 ^ **• 


^” B nn f n^ 

< B irA> 


[t } = 


(B- 2 e 2 ) 

0 


[Ml 


ini =: 


0 

0 


• ^ B nm e m^ 


0 


0 

0 


0 


(B n,m+1 e m+1* # * 

(B m+1,m+1 e m+1^*' 


^ B n,m+1 e m+1^ 




(B" -.f . ) ... 

v n,m+1 m+1 


< B An e n> 


K^,n e n y [ 


(B" e ) 
v nn n' 


^ B m+1 ,n f n^ 


(B 51 f ) 
' an n' 


2e, 


0 


0 


0 


[ S] =* 


0 


2 e. 


m 


0 


0 


IT] 


2fr 


0 


0 


2 f 


m 


0 


0 


0 


0 
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In the above derivation, the slack bus is numbered as 1, 
P-V buses as 2 to m and P-Q buses as m+1 to n. To achieve 
decoupling, the following further assumptions are made : 

i) A flat voltage start (e+jf = I.O+ 3 O.O) is used. 

ii) The matrices obtained using the flat voltage start are 
treated as constant at every subsequent iterative step. 
This avoids the re-calculation and re-tri angulation of 
the Jacobian matrix at each Iteration. The algorithm 
obtained, therefore, utilizes the fixed tangent slopes 
corresponding to the flat voltage start instead of the 
varying tangent slopes at each iteration. 


With the above assumptions, equation (2.22) is reduced 
to the following set of decoupled equations in rectangular 
coordinates : 



<0 


* 




[ B’l 





[b !, i ■ . ; 


* > 



[ B”' ] 





(2.23) 


(2.24) 


(2.25) 
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■where [ B’"J = 2 f I | 

Since [ Af 1 and hence [f 1 are known for the P-Y buses 
from (2.23), [el for the P-7 buses can be calculated using 


its specified voltage magnitude from 

[e 2 1 = [V 2 J - [f 2 1 (2.26) 

Replacing (2.25) by (2.26), the final PDLR equations become, 
[ AP r J = [ B' 1 [ Af ] (2.27) 

(el = (IV 2 ! - Tf 2 !)^ (2.2S) 

[ A Q r l = [3"1 [ A e 1 (2.29) 


The above equations can also be obtained from RM load flow 
equations in rectangular coordinates as indicated below. 

2.4.2 Derivation of PD1R from RM load flow equations 

Using the flat voltage start tangent slopes at every 
iteration, the elements of the Jacobian in RM equations 

(2.11) [Ji 1 (vide Appendix A) are modified as 

l— ' j ^ j ••• n* 

follows : 


3 P 


r 

i 





3pF 3 qF 

1 _ _ -n 

3 -B iq 


(2.30) 


(2.31) 
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»* 

: r 'jrn » «> » 

aC j 


a i + °ii 


* 

“if. 


\ - B ii 


3 e* 


-b. - 33* - 

l ii 


3 Q? 

l 

a fV 

i 


a. - G. . 

l ii 


where 


a. = 


n 

3=1 10 


n 


\ 


l B 

i=1 


( 2 . 32 ) 

(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 


Prom equations (2.36 and 2.37), we may assume a-^ << G.^ and 
i <<: hence a^ and b^ can be neglected in equations 

(2.32 - 2.35). Further, neglecting G^ and G.^ in the above 
expressions, equation (2.11) is reduced to 


[AP r ] = 

'V 

[Af ] 

t AQ r l = 

[B Q ] 

[Ae] 

[Ay 2 ] = 

t B"'] 

[ Ael 


(2.38) 

(2.39) 


(2.40) 
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•where (B^I and [ Bq J are the elements of [-B 1 matrix. 
Changing the matrices [B ] and [Bq| to [B’J and [B" ]with 
the help of the modifications suggested in reference £291 and 
replacing (2.40) by (2.26), we get the final load flow equa- 
tions in rectangular coordinates which are the same as obtained 
from BD1P by linear transformation in equations ( 2.27 - 2 . 29 ). 

2.4.3 Iterative scheme 

Equations (2.27 - 2.29) are solved successively using the 
most recent values. The best iteration strategy is termed as 
the (If, 1e) scheme, which comprises one solution for [ Af J 
using equation (2.27) to update [ f ] relating to the P-V 
and P-Q buses, one solution for f.e j using equation (2.2s) 
relating to the P-V buses and then one solution for[Ae] using 
equation ( 2 . 29 ) to update [ e} relating to the P-Q buses. 
Iterations are repeated until the real and reactive power mis- 
matches are within a specified convergence tolerance. In the 
solution of lAf] andfAeJ, the constant sparse triangular 
factors of tB’ } and [ B" ) are used. Interchanging equation 
(2.28) with equation ( 2 . 29 ) results in poor convergence. 

2.4.4 Consideration of Q-limits of voltage controlled buses 

To take into account the Q-limits of the voltage contro- 
lled buses, the following two methods can be used. 

i) Method of switching P-V bus to P-Q bus and vice versa. 
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Each violated P-V bus is explicitly converted to a P-Q 
bus to fix the reactive power at the limiting value. The com- 
putations are done as for a P-Q bus and this requires the re- 
triangulation of [ B n j * If in the subsequent computation, the 
reactive power falls within the available value, the bus is 
switched back to a P-V bus at the specified voltage. 

ii) Method of sensitivity factor 

The change in e^ at the ith violated P-V bus is calculated 
by using the sensitivity factor [29] given by 

Ae i = S i AQ i 

This is used to reduce the error to zero. 

The sensitivity factor may be computed by augmenting the 
existing factors of [B u ] as indicated in Appendix B. This pro- 
cedure does not require the re-tri angulation of [B"] . 

2.5 LOAD PLOW SOLUTION BY NEWT ON-RI CHARD SON METHOD (HEM) 

2.5.1 Derivation of NKM load flow equations 

The basic idea of this method is to combine both the Chord 
method (or simplified NM [12])and the Newton iterative method 
[131 . The linear set of equations appearing in NM has been 
represented in a matrix form in Section 2.3 (vide equation 2.6) 
as 


ted* 1 > [yi )i 


[ j 1 [AX ] 


(2.41) 
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Let the Jacobian [ J ] at the kth step or major iteration be 
split as 


rn (k) 1 

[J1 =■-■{£) 


tj) ^ 


(2.42) 


where 


[J] = Jacobian evaluated using the initial bus voltage 

magnitude and phase angle estimates 


[j/ k) = I J] (0) - tJJ (k) 

a 


a = real constant 


Using this splitting, the Gauss-Seidel iteration for solving 
equation (2.41) at the kth major iteration is written as 

' Tk) 1 Jl (0) tAxl (k) = t J 1 ] fAxj (k) + [g([x], [y] ) ] (k) 

a i + 1 i rt 

0 0 (2.43) 

where the subscript i Q represents the minor iteration. 
Combining (2.42) and (2.43) and simplifying, we get 

™ Xk) 1 J l(0) U d (k) = “ I' Jl (k) l Ax]^ k) + [ g([x] , [y] )} (k) 
a i o 0 

(2.44) 


l A xj ^ =(4x]^ + [Sx]^ 

i +1 i i 

0 0 o 


(2.45) 
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^ ir ^ 

A sequence of minor iterates [ fix]^ J > CAxJ-j , [Ax]^ ••• 

[ Ax]^ k ^ is generated using equations (2.44) and (2.45) to 

approximate the solution of equation (2.41) at the kth major 

fk') 

iteration in which [Ax]^ =0. 

(k) 

The final converged value of [ A xj v at the kth major 
step is given as 

[Ax] ^ = [5 x]^ + [<5x}p) + [Jx]^ + ... (2.46) 

where 

K 

[«x]£ k) = a (k) [J] (0) {[g([x] , [y]) ] (k) } 

[6x]{ k) = a (k) [J] (0) { [g([x] ,[y] )] (k) - [ J] (k) [ 6 x]^ k) > 

[ 6x]^ k) = a (k) lJ] (0) { [ g([x] ,[y] )] (k ^ - [J]^[[ 6x]^+{6x3^] 
and so on. 

A choice value of a between 0.0 and 2.0 is indicated in [37]. 

We have found that a value a = 1.0 gives good results for the 
different cases studied. 

2.5.2 Iterative schemes 
(k) 

let m^ ' be the number of terms included in equation 

(2.46) at the kth major iteration. Choosing m Q = 1 and a =1 

at each major iteration reduces the proposed method to the Chord 

or simplified MM method. The numerical results based on three 

(k) 

iterative schemes depending on the value of m^ at each major 
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iteration have "been considered in Section 2.7. Table 2,2 
gives the pertinent information regarding the three schemes. 

TABLE 2.2 ITERATIVE SCHEMES OF HRM 
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The "bus admittance matrix [Y] is assembled in such a manner 
that the first n^ buses represent the voltage controlled 
buses including the slack bus and the last n^ buses repre- 
sent the load buses. 

The equivalent load admittance at bus i is given by 



P S P _ i oft' 
±ji 3 ni 


Sp 



(2.48) 


If at each load bus, the scheduled complex load power is - 
represented by the equivalent admittances, then the vector 
equation (2.47) gets modified to 


t I 1 = [Y’l [E ] 


(2.49) 


i.e. 


j 

i. — 


toj 


; t t g G i 


I :i LG 1 


1 y gl> : ;■ [ i 


I y LLl ! ! t%]j 


<- ! 

2 
V 

X 

9 

I 


and 


£ Y ;W - + * 


(2.49a) 


where 

1 Y LL 1 = load bus admittance matrix before adding the 

equivalent load admittances. 

[ ] = diagonal matrix representing the equivalent load 

admittances. 
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From (2.49), we get 

[isl - {tr GG l - tr sl ] ■ 1 t^ G l ) 

[EjI = - tY^r 1 tr LG l [E 5 1 

Equations (2.50) and (2.51) can be written as 


(2.50) 

(2.51) 


i - <*4s> tE &> (2 - 52) 

[Y iL l (Ej-1 = - [7 lg 1 [E q ] (2.53) 

where [Y^q.] represents the bus admittance matrix of an equi- 
valent system consisting of only voltage controlled buses and 
the slack bus. CYq.q.1 is obtained by Kron's reduction process, 
which' essentially consists in applying Gaussian elimination 
to the modified bus admittance matrix tY' 1 , exploiting the 
sparsity of the network as explained in [21,24], At the end 
of the reduction process, the form of the modified bus ' 
admittance matrix [Y 1 ] is shown in the accompanying Fig. 2.1. 
The reduction process provides the triangular factors of 
[Y| Ti ] matrix in addition to the non-triangularized equivalent 
admittance matrix [Y^q.] as indicated in Fig. 2.1. 

let 

CT&jI . = t g r) + 3 t B r ] (2.54) 


The well known polar power mismatch NM is applied to an equi- 
valent system consisting of only voltage controlled buses. Efy 
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approximating the resulting Jacobian with the help of the 


valid assumptions [ 29 ] namely, cos 6. . « 1 and 
G Rlj Sln 6 ij << B Rij’ we get 


‘ V 


- ^ U 8 v ) 


(2.55) 


where the elements of [B’ ] are strictly the elements of 

t- B r ] 


ip . 
vi 


real power mismatch at the bus i relating to an 
equivalent system consisting of only voltage con- 
trolled buses and the slack bus 


n 


G 


- I {V 3 (G Rlj 003 9 i3 +B Ri5 sln 9 lj )) 


3=1 


r R ij 


and 


[ A 8 ] 

v 


voltage phase angle correction vector relating 
to the voltage controlled buses. 


The final load flow equations are represented by equa- 
tions (2.55) and (2.53). The technique to solve for 
from equation (2.55) is by KRM (vide Section 2.5) without 
t ri an gul arizing [B^J at every iteration. 

let [B^ o J be the matrix obtained during the first 
iteration. Then, the final converged value for [A8 V J at each 
iteration k is expressed in a series form as 
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[ie v ] (k) = t 66 vl' k) + tse j/ k) + t«e jW + .:<2.56) 

c vi v 2 


■where 


t6e v } (k) = IBy _1 { t 4 ^-i } 


t «e l (k) = (B’j r 1 {[ A ^ v ] - [B^j (k) t«e v ! (k) } 

' V o 


and so on. 

The experimental studies indicate that it is only necessary 
to include the first two terms of the series given by (2.56) 
i.e. ITEM Scheme I (vide Section 2.5) for the same solution 
accuracy, as obtained by the direct triangularization of [B^J . 
Hence, each solution of C AQ^.1 from equation (2.56) requires 
two repeat solutions using the factors of [B-j^ ]. 

2.6.2 Iterative algorithm 

The basic iterative algorithm for the proposed method is 
as follows : 

1) Read the system data and initial estimates of the bus 
voltages (magnitude and angle). 

2) Borm the bus admittance matrix [Y ] using the network 
elements and shunt elements except loads. 

3) Obtain the equivalent load admittances corresponding to 
the scheduled complex load powers using equation (2.48) 
and modify [Y J to get [Y'] . 
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4) Obtain the equivalent bus admittance matrix [YqqJ by 
reduction process and hence [B^]. 

5) Obtain the correction [A8 V ] using MM Scheme I 

. (equation 2.56) and determine the new values of the 
complex voltages relating to the voltage controlled 
buses. 

6) Obtain the load bus voltages using the triangular 
factors of [Yj-^] available from step (4) in equation 
(2.53). 

7) Steps (3) to (6) represent one complete iteration. 
Iterations are repeated until the exact solution is 
obtained, i.e. when 

- ' Max | AP vi | < e t 

l 

The speed of computation can be further increased by 
a modification described below. 

2.6.3 Modified form of NETR 

The above algorithm requires more computing time than 
EDIiP due to the necessity of reducing the modified bus 
admittance matrix [Y* ] at each iteration, even though, the 
total number of iterations are less for moderately accurate 
solutions (mismatch of 1 MV/). To reduce the computation 
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time significantly, in every alternate iteration, the 
changes in [D-j-J matrix in equation (2.49a) due to the 
changes in load hus voltages are neglected for the calcula- 
tion of However, the changes in [D^] matrix for 

the calculation of the load hus voltages from equation 
(2.53) can he incorporated by writing it as follows. 

t r il> IE l’ IE L 1 - W 1G 1 (Eg! (2.57) 

where the change in iHj] is represented by 


t ABjl 


old new 


(2.58) 


This allows the old triangular factors of [ Y { Ti } to be used in 
the calculation of 

2.6.4 Consideration of Q-limits of voltage controlled buses 

Each violated voltage controlled bus is explicitly con- 
verted into a load bus, so that the Q-limit is held at its 
specified value. If a P-7 bus is switched to a P-Q bus, then 
an extra row and column is eliminated from HqqJ relating to 
the switched bus in order to obtain a new equivalent system 
consisting of only non-violated voltage controlled buses. The 
complete reduction of [ Y' ] matrix is not necessary once again 
at any iteration due to the violation of voltage controlled 
buses. If in the subsequent iteration, the reactive power 
falls, within the specified limit, the bus is switched back as 
a P-7 bus at the specified voltage. 



51 


2.7 NUMERICAL RESULTS ON TEST SYSTEMS AS 33) COMPARISON 

The proposed load flow methods outlined in Sections 
2.4 to 2.6 have been tested by computing the load flows for 
the following power system examples on DEC 1090 system. 

1. IEEE 14 bus test system 

2. IEEE 14 bus test system with R/X ratio of the heavily 

loaded line increased by 5 times 

3. IEEE 30 bus test system 

4. IEEE 57 bus test system 

5. 28 bus heavily loaded system [38 ] 

The specifications of the above systems are given in Appendix D. 

In addition, EDLR has also been tested on a 133 bus system given 

in reference [22], 


The results obtained from the load flow solutions on these 
systems using the flat voltage start are listed in Tables 2.3 
to 2.10. Details regarding the max f AP | and max | A Q | f or 
each iteration are given in Tables 2.3 to 2.6. The number of 
iterations and the relative CPU time in 1 ° are shown in Tables 
2.7 to 2.9. Tables 2.3 to 2.9 pertain to unadjusted solutions 
(no Q-limits imposed) so as to compare the convergence charac- 
teristics of the proposed load flow techniques, whereas 
Table 2.10 corresponds to adjusted solutions with the considera- 
tion of Q-limits. Eor the consideration of the Q-limit vio- 
lations in EDLP and EDLR, method 1 , discussed in Section 2.4.4 


is used. 


Eor comparison, the load flow^esults obt.p|iu| 


&£« ft*. 


■m ^ 



TABLE 2.4 MISMATCHES DURING DM, NRM, EDLP, EDLR AND NETR LOAD FLOW SOLUTIONS 
IEEE 30 BUS SYSTEM 

MISMATCH < 0.01 P.U. MW AND 0.01 P.U. MVAR 
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TABLE 2,5 MISMATCHES LURING NM, NRM, FDLP, PDLR AND NETR LOAD FLOW SOLUTIONS 
IEEE 57 BUS SYSTEM 

MISMATCH < 0.01 P.U. MW AND 0.01 P.U. MVAR 
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TABLE 2.6 MISMATCHES DURING- IDLE AND PDLP LOAD PLOW 
SOLUTION 

133 BUS SYSTEM [ 22 ] 

MISMATCH < 0.01 P.U. MW AND 0.01 P.U. MY AR 


■Iteration * Max J A P | 
!No. i. __ 

' Hedlp" 

t f. 


; 0 I 4.2352 

s 1 ; 2.2161 

i 

] 

j 2 ■ 0.1060 

* i 

■ 3 I 0.0070 


mismatch in p.u. ; Max |AQ {mismatch in 



p.u. 


“"PDLR ’ r ~ 

! 

PDLP ' ; ” 

! 

PDLR ■ 

4.2852 : 

26.3119 ' 

26.2690 \ 

2.0897 'i 

0.1732 j 

0.1385 , 

0.0445 j 

0.0085 ; 

0.0127 l 

0.0022 


0.0020 ; 


TABLE 2.8 COMPUT AT IONAL COMPARISON BETWEEN PDLR AND PDLP 
133 BUS SYSTEM [22] 


"Method 


| PDLP 
JPDLR 
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MISMATCH < 0.01 P.U. MW AND 0.01 P.U. MY AR 



No. of 
iterations 



Relative computer j 

time in $ on DEC 1090 ] 

syst em i 
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TABLE 2,7a COMPUTATIONAL COMPARISON BETWEEN PELP, EDLR AND NETR LOAD PLOW SOLUTIONS 
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T ABLE 2.9 OQMPUT AT ION AL COMPARISON B2FWEEJ? NM AUD NRM LOAD PLOW SOLUTIONS 
MISMATCH < O.OOl P.U. MW AND O.OOl P.U. MVAR 
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TABLE 2.10 CONS IDEM ION OE Q-LIMITS 
IEEE 30 BUS SYSTEM 




MISMATCH ■< 0.1 MU 

AND 0.1 MVAR 

/■v— ■ > .1 

[Method 

l 

1 

l 

b 

No. of iterations 

* 

Relative computer; 
time in °f° on j 

DEC 1090 system j 

i 

i 

[edlp 

1 

l 

l 

4* ; 

O 

O 

sEDLR 

i 

k 

1 

64 1 

» 

105 J 

t 

>NETR 

I 

5 

l 

4 ! 

117 I 


f 

» 



NM and EDLP are also included. The same ordering and pro- 
gramming techniques are used for all these methods. The 

(]r) 

value of o v = 1 is used for NEM load flow solution. 

Based on the results shown in Tables 2.3 to 2.10, the 
comments and comparisons made with respect to the proposed 
load flow techniques are given below. 

A. Comments on EDLR and its comparison with EDLP 

1) The simulation time required per iteration of EDLR is 
smaller than the corresponding iteration of ■ EDLP. 

2) Although EDLR takes more -number of iterations for con- 
vergence as compared to EDLP, the total computer time' is 
less by about 5 to 10 percent for all the test systems 
(no Q-limits imposed) 



60 


3) For the IEEE 14 bus system with E/2 ratio of the 
heavily loaded line increased by 5 times, the saving in 
computer time by FD1R is more than 20 percent. 

4) ^or the heavily loaded 23 bus system given in reference 

[38] , with some of the phase angles at the system 
buses of the order of -55°, both FDLR and FDLP took 
considerably more number of iterations for convergence 
and the saving in computer time by FDLR is about 13$. 

5) With 2 buses violating their Q-limits on the IEEE 30 bus 
system, EDLR took 2 iterations more than EDEP. However, 
the total computer time is comparable with FDLP. 

B. Comments on HEM and its comparison with EM 

1) The rate of convergence is adjustable and depends on 

k 

the value of m . 

o 

2) The solution accuracy and the convergence rate of Scheme 
II are similar to those of EM. 

3) The computational time required by all the schemes is 
less than by EM and it decreases, as the size of the 
system increases. For the IEEE 57 bus system, the com- 
putational time saving is 36 percent. 

4) It is known that in EM, the convergence rate is not 
affected when R/X ratio of the transmission line is 
increased. This was also found to be true in HEM, when 
R/X ratio of the heavily loaded line was increased by 

5 times in the IEEE 14 bus system. 
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5) Eor the heavily loaded 28 bus system, with some of the 

bus angles at the system buses of the order of -55°, it is 
essential to use 3HM after two initial NM iterations to 
obtain successful convergence. 

C. Comments on ESTR and its comparison with EDLP 

1) The solutions are sufficiently accurate in 2 or 3 iterations 
for all the test systems. The number of iterations is not 
seen to be a function of the total number of buses. The 
total computing time is more as compared to EDLP , but is 
much less as compared to NM. 

2) The convergence rate is not affected on the IEEE 14 bus 
system, with R/'l ratio of the heavily loaded line increased 
by 5 times and also on a heavily loaded 28 bus system. 

3) The modified form of the algorithm which ignores the changes 
in load bus voltages in alternate iterations in the calcula- 
tion of [ A Q v } is efficient, faster and requires comparable 
or less computing times compared to EDIP. 

4) The total simulation time is comparable to EDIP, on the IEEE 
30 bus system with two buses violating their Q-limits. 

2.8 COECIUSIOIT 

In this chapter, three new load flow solution methods have 

been presented, viz. 1) fast decoupled method in rectangular 

coordinates possessing distinct computational advantages over 
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the polar version, 2) Eewton-Richardson method, an improved 
version of JJM and 3) network reduction method, a comparable 
alternative to the fast decoupled method for medium size 
systems. The IDLE, and ITETR' can be easily extended for rapid 
contingency evaluation needed in power system security assess- 
ment as shown in the next chapter. Besides, 3TRM and NETR load 
flow algorithms can be efficiently adapted respectively for 
optimal load flow and economic dispatch calculations as illu- 
strated in Chapter IV. These methods have been tested on the 
IEEE 14 bus, 30 bus, 57 bus test systems plus a heavily loaded 
28 bus system and a 133 bus system. 



CHAPTER III 


CONTINGENCY EVALUATION 

3.1 IHTRODU CT ION 

In Chapter I, we presented the basic security concepts 
and the various security related functions. The key function 
of the security analysis is known as contingency evaluation 
which identifies, whether the normal system is secure - 
or insecure. A number of fast on-line approximate algorithms 
[20,57-59] have been proposed for contingency evaluation. 

The work reported in this area along with a critical review 
of the pertinent literature has been presented in Chapter I 
(Section 1.5.2), The techniques used for contingency 
evaluation must possess the essential features of fast 
computational speed, reliability, reasonable accuracy and 
moderate memory storage requirements. The objective of this 
chapter is to extend EETR & PDLR load flow techniques 
developed in the preceding chapter for contingency evalua- 
tion. The proposed techniques make use of only base case 
factors for the simulation of different contingencies. 
Detailed solution algorithms using the proposed techniques 
are given for implementation purposes. The proposed algo- 
rithms have been verified by simulation on the IEEE 14 bus, 

30 bus and 57 bus test systems. A comparision of the pro- 
posed methods, and a competitive method, EDLP is made in 
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terms of solution accuracy, core storage and computational 
requirements. 

3.2 STATEMENT OP THE PROBLEM 

Contingency evaluation consists of two parts. Firstly, 
it determines the steady state solution of the nonlinear load 
flow equations due to the assumed network and/or power outages 
as classified in Table 3.1. Secondly, it analyses the results 
for different constraint violations like transmission line and 
transformer power flows, reactive power and bus voltage limits. 
The thesis is mainly concerned with the first part of the pro- 
blem, which can be restated as follows : 

Given : the base case load flow solution of the system 
along with the contingency list. 

To determine : the post contingency states of the system. 

TABLE 3.1 CLASSIFICATION OF OUTAGES 


\ Type of outage 

J 

* Network 
\ 

$ 


X 

1 

s 

> 

\ 

i 

1 


Power 




; ii) 

3 

j in) 

i 

i 



£ . . . \ 
> 111 ) 

i 


System status \ 

\ 

' " * "1 

4 

outage of a transmission line J 

, j 

outage of a transformer * 

outage of a shunt capacitor ? 

or reactor 

loss of a generator 

sudden loss of a load 

sudden change in flow in 
inter -tie 
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5.2.1 Contingency list 

Since it is very difficult to test a large power system 
for every conceivable disturbance, a suitable contingency 
list is generally created based on the statistical analysis 
of past occurrences, off-line simulation studies etc. The 
contingency list can be prepared in a systematic way to con- 
tain only the cases pertaining to the present state of the 
system. It should be updated frequently according to the new 
information available. More recently, Ejebe and Wollenberg 
[62J have proposed a fast technique for automatic ranking and 
selection of contingency cases. The algorithm developed, 
ranks the line and generator outage cases according to the 
severity of their effects oh bus voltages or line flows using 
the d.c. load flow as -well as nonlinear a.c. load flow equa- 
t ions . 

3.3 CONTINGENCY EV-ALU.A1I0N BY NETWORK REDUCTION MW NEWTON- 

RJCHARDSON TECHNIQUES 

3.3.1 Development of solution technique 

Eor the evaluation of network outage between the voltage 
controlled buses including the slack bus termed as type-1 , 
only 4 elements of the base case [Y^J and hence of and 

[B^ ] (vide Chapter II, Section 2.6) are changed. If the outage 
occurs between the buses i and then the change in the re- 
duced matrix JY^qJ is 
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< ay gg-1 


i ! -y _• j - 


I -y 

i 
> 

I 


ID 


10 


ID 


ID 


I “ Y iD 


- y 


DO 


(3.1) 


■where 


y. , y. = line charging capacitance and shunt branches of the 

-L w J U 

off-nominal transformer equivalent circuit lum- 
ped at nodes i and 3 

Hence, [Y^q. ] and [ ] are modified as. 




new 


- '%) + 


[ B^l 


new 


T 


t B Rl + ‘,IV f«v 1 


where 


and 


= b . . 
ID 


[ } is a highly sparse vector given by 

1 . D 

t0,0, ..., 1,0, ... -1,0, ... 0] 


(3.2) 

(3.3) 


(3.4) 


If the network outage involves the removal of an off-nominal 

transformer, then vector becomes 

3 ; 

[0,0, ... i , 0, ... -1,0, ... 0] (3.5) 

iD ! 
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By using the inversion of a modified matrix method (vide 
Appendix B), the correction for[A9 v ] is obtained as, 

- V V "V 1 46 ^ 1 (3,6) 

where 

-1 m 

txy = tsy uy 1 (3.7) 

and 

a, - (£- + [jgny r 1 (3.3) 

Bor each network outage of type 1, the non-sparse vector ty 
is determined using the factors of [By in equation (3.7). 
After each solution of [A9 V J using equation (2.55), it is 
corrected by an amount using equation (3.6). After one or 
two iterations, depending on the needed accuracy, the load 
bus voltages are determined using equation (2.53). 

When the network outage occurs between a load bus and 
any other bus or between the load buses termed as type-2, 

[Yjy itself changes in equation (2.49). Without going 
through the reduction process once again, [Y^q. 1 and hence 
[By is obtained using the matrix updating algorithm pro- 
posed by Smith [21J (vide Appendix B) which makes use of the 
base case factors only. 

[ in equation (2.55) is solved by WBM-Scheme I 

(vide Chapter II, Section 2.5) with two repeat solutions 
using the factors of [BL ) relating to the base case. Further, 
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since f Y ' J also changes at the maximum by 4 elements, the 

JjJj 

lo ad bus voltages are determined using the inversion of a 
modified matrix' method. 

After the solution of (2.53), [ Ej is corrected by an 
amount 

= - a^l X L ] (ty ty (3.9) 

where 

ti y T (3 - 10) 

y = (jj + ny ly ) _1 

and 

^1 = “ y iD 

Details of the derivation of vectors [ and [ and scalars 
and are given in the Appendix B. 

3.3.2 Impl ement at i on 

The technique described in the preceding section for 
contingency evaluation is now given in the form of a detailed 
algorithm for implementation purposes. 

Iteration scheme 

Two iteration strategies with different speed and accu- 
racy are proposed. In scheme I, termed as (10 v , scheme, 

first 0 V - iteration is performed using equation (2.55) 
followed by 1^ - iteration using equation (2.53). This is 
found to work well for most of the outages. However, for the 
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outage of a very heavily loaded line, scheme II termed as 
(2e^. f lE^) scheme is proposed. In this scheme, which gives 
"better accuracy, 2 q ^ - iterations are performed first and 
then followed "by one - iteration. 

Preliminary step 

We perform load flow solution for the base case and 
store the triangular factors of [Bi] , [ Y , the factors 
needed for the matrix updating algorithm, the matrix fY GG ] 
and the converged bus voltage vector. Por on-line studies, 
the bus voltage vector may be obtained from the state esti- 
mator and the different matrices may be easily built and 
factorized. At the end of the base case load flow solution, 
the table of factors and the various matrices are stored as 
shown in the accompanying Pig. 3.1. 

The main steps involved in the iterative solution are 
listed below : 

A. Por a network outage between voltage controlled buses 

including the slack bus - type 1. 

1. Obtain [ using equation (3.1) and form [ Y 

A? v new 

2. Compute the mismatch [ 

v 

3. Computet A0^] by a repeat solution using the base 
case factors of EB^] in equation (2,55) 

4. Compute the correction for t using equation (3.6) 




FIG- 3.1 STORAGE OF DIFFERENT MATRICES AT THE END OF BASE 
CASE LOAD FLOW SOLUTION • 
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5. For (10 V » lE^) scheme, proceed to the next step. 

For (29 v , IEjO scheme compute the mismatch using the 
new values of [ ] and repeat steps (3) and (4) and 

proceed to step (6). 

6. Compute {- [ [Eg] } and solve for using 

equation (2.53) 

7. Repeat steps (1) to (6) for each network outage of 
the above type to be simulated, 

B. For a network outage between a load bus and any other 
bus or between load buses - Type - 2. 

1 . Obt ain [ A Y GG- J using the matrix updating algorithm 

and form [Y^q.] and hence [B^] 

new new 

aY v 

2. Compute the mismatch [ -y ] 

v v 

3. Compute [A©^] by NRM - scheme - I by two repeat 
solutions using the base case factors of [ B^] and 
update the angles [0 V J 

4. For (10 V > lE^) scheme, proceed to the next step. For 
(20^, lE^) scheme, compute tho mismatch using the new 
values of I 0 V J and repeat auiStcp (3) and proceed 

to step (5). 

5 . Comput e {Y LG } [Eg] } and solve for py using 
equation (2.53) 

6. Obtain the correction for { using equation (3.9) and 
update the voltage vector { E^l obtained from step (5). 
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7. Repeat steps (1) to (6) for each network: outage of the 
above type to be simulated. 

G. Power outage 

1. Redistribute the power lost according to the generator 
distribution factors and proceed to step (7). If it 
is not possible to do it, proceed to the next step. 

2. The real power generation lost is automatically dis- 
tributed to all the other regulating generators on the 
basis of their capacities and governor droops based on 

the steady-state frequency drop of the power system. 

AP v 

3. Compute the new mismatch [ } 

V 

4. Compute [ A8 v ] using equation (2.55) and update the 
angles { 0 y } 

5. For (19 V , lE^) scheme, proceed to tha'noxt _ step. . For 
(^ v , lE^) scheme, repeat steps (3) and (4) and 
continue. 

6. Compute {- [ YjJ [ E & ] > and solve for [Ey using equa- 
tion (2.53) 

7. Repeat steps (1) to (7) for each power outage to be 
simulated. 

3.3.3 Solution speed and storage 

A realistic comparison is made in Tables 3.2 and 3.3 for 
the major computing operations to be performed for the simu- 
lation of outages using the proposed method and the FDIiP, 
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TABLE 3.2 MAJOR COMPUTING OPERATIONS NEEDED EOR 
6 (EDLP) AND ! 6 V (NETR) - ITERATION 


EDLP 


Proposed- NETR 


(Operation 

\ 

i ( 1 e , iv) 

(20,2V) 

i 

i 

f 

Oe v , i Bj,) 

(2e v , IE 

i 

I 

^ v ur - * 

• 

i 

, ... 


fBase case matrix 

j IB* 1 

[B'l 


lB k J 

[to be built and 

i 


f 

.a. 

[factorized 

y 

;(n s xn g ) 

(a s m J 

5 

ft 

K xn v ) 

(n^xr^) 

\ Matrix to be up- 

| 

— 

t 

i 

« b r1 

pb) 

-dated per type-2 
l out age only 

| 





« 



< 

\ Number of mis- 

i i 

2 

£ 

1 

1 

2 

match computa- 

> 


t 



tions per outage 

$ 


a 

i 



‘Number of repeat 
• solutions per 

\ 2 

3 

J 

\ 

2 

3 or 4 

> 





jout age 

i 

t 


f 




The following comments are made from the Tables 3.2 and 

3.3. 

The schemes (I 0 , IV) and (l0 y , 1E^) require the same 
number of mismatch computations and repeat solutions. The 
(20^, lE^) scheme needs one mismatch computation and one 
repeat solution less than the (20, 27) scheme. The proposed 
schemes require computing operations for updating the matrix f 
[ for every type-2 outage. The arithmetic calculations 
needed for { -{lpQ.1 [E Q ]} computation relating to the proposed 
schemes are less than the corresponding mismatch computation 



74 


TABLE 3.3 MAJOR COMPUTING OPERATIONS NEEDED EOR V 
(EDIP) AND E l (NEER) ITERATION 


[Operation 


< EDLP 

I (1 9, IV) (26, 2V) 


[Base case matrix : [ B ,! J 

Ito be built and i 

ffactorized l (n^xn-j-) 

l i 

l Number of mis- j 

j match computa- ; 1 

*tions per outage { 

: i 

> Number of | 

.computations per [ 

\ out age ■ 

' f % 

\ i 

i Number of repeat \ 

^solutions per j 2 

•out age j 

\ i 


[ B" ] 

(n^) 


2 



; 1 

! Proposed -_NETR l 

J Ov 1E I> < 2 V 1E L* 

* 5 

t " ' ■ ' ‘ : ‘ ‘ 

| t Y ii) t r id I 

j (njmj) (^nj) i 



requirements of the EDIP schemes. Further, the total computing 
time needed for - iteration relating to the proposed schemes 
depends on the number of P-V buses in the system, since 8 - ' 

iteration is performed on the reduced system consi sting of only 
P-V buses. 

A numerical comparison reveals : 

i) Eor type-2 outage, the proposed method is faster than EDIP, 
if the number of P-V buses is less than about 10$ of the 
total number of buses in the system. 
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ii) For type-1 outage, the proposed method does not require 
the updating of [B^] . Hence, the reduction in time as 
given in (i) is increased hy about two times. 

The matrices involved in the proposed method are symmetric 
and hence, only the upper or lower triangular factors need to 
be stored for carrying out repeat solutions similar to FDIiP. 

3.4 CONTINGENCY EVIL UATION BY FAST DECOUPLED TECHNIQUE IN 
RECTANGULAR COORDINATES 

3.4.1 Development of solution technique 

With the help of the inversion of a modified matrix 
method, (vide Appendix B) network outages like lines or trans- 
formers are simulated. This method avoids the re-tri angulation 
of the network matrices [B 1 ] and [B”J for the simulation of 
single or multiple outage cases. Experimental studies indicate 
that it is only necessaiy to simulate the removal of series 
transmission elements from these matrices. However, all the 
changes must be correctly taken into account in the calculation 
of active and reactive power mismatch vectors. 

For each network outage involving lines or transformers, 
two vectors [X 1 ] and tX”J are calculated in the following 
manner. 


I X* 1 

= ( 3*1 

T 

CM’] 

(3.12) 

tx !i j 

= IB' 5 ]" 1 

T 

[M*1. 

(3.13) 
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where [Ml and [M”l are highly sparse row vectors each 
containing one or two non-zero elements (vide equations 3.4 
and 3.5). 

After each solution of equation (2.27), {A fl is com- 
pensated "by an amount 

- a' [ I’llM'] [Af ] (3.14) 

Similarly, after each solution of equation (2.29), [AeJ is 
compensated by an amount 

- a" t X”] [M ! >] [Ae } (3.15) 

where 

a' = (p + [ M’J [X>] ) -1 
a” = (p, + [ M<‘] l-Xl )~ 1 

and 

b’ , b" = susceptances of the outaged line or transformer. 

There is no correction needed for the equation (2.28). De- 
tails of the derivation of vectors [X'] and [ X" 1 and scalars 
a* and a” are given in Appendix B. 

3.4.2 Iterative algorithm for outage simulation 

A. Base case load flow 

1. Read the input system data and initial estimates for 
[e) and [ f ] 

4 

2. Form the bus admittance matrix [Y] 
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3. Obtain IB*] and [B i! ] from [ Y 1 using the approximations 
and modifications suggested in reference [291 (vide 
Chapter II, Section 2.3.2) and store its table of 
factors 

4. Obtain the base case solutions { e Q J and f f 0 3 "by 
iterating the equations ( 2.27 - 2 . 29 ) successively 
until the convergence is reached and store it. 

B. Hetwork outage 

1. Initiate the base case solution 

2. form the vectors [M ! J and [M" ] corresponding to the 
outage case and solve the equations (3.12 and 3.13) to 
obtain [X’J and [X” ] 

3. Compute [AP r J using equation (2,12) and solve for [Af J 
by a repeat solution using equation (2,27). 

4. Compute the correction for [Af ] using equation (3.14) 
and update [f I 

5. Compute [e } relating to the P-Y buses using equation 

( 2 . 28 ) 

6. Compute [AQ r ) using equation (2.13) and solve for (Ae J 
relating to the P-Q buses by a repeat solution using 
equation ( 2 . 29 ) 

7. Check whether the network outage is between P-Y buses 
including the slack bus. If it is, proceed to step 
(9). Otherwise, continue to the next step. 

8. Compute the correction for [AeJ using equation (3.15) 
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9. Update ( e ] 

10. For (If, 1e) scheme, proceed to the next step. 

For (2f, 2e) scheme, repeat steps (3) to (9) and 
continue. 

11. Repeat steps (1 to 10) for each network outage to he 
simulat ed. 

0. Power outage 

The power outage can he simulated based on thelMmilAr lines 
described in Section (3.3. 2-0). 

3.5 HUMERI CAL BESTED) S OR TEST SYSTEMS MB COMPARISON 

The IEEE 14 bus, 30 bus and 57 bus test systems (system 
specifications are given in Appendix D) were used for contin- 
gency evaluation with the proposed NETR and FDLR techniques 
and the details of the results are compiled in Tables 3.4 and 
3.5. The number of line/transformer outages simulated in the 
three test systems were 19,26 and 35 respectively using 
(10 Y , lE^), (20 y , lE^), (If, 1e) and (2f, 2e) schemes of the 
proposed methods and (1®, IV) and (20, 2V) schemes of FDLR. 

All the line and transformer outages were considered for the 
14 bus test system, excepting line 10, whose outage splits the 
system into two networks. Most of the lines selected for the 
other systems were heavily loaded. Table 3.4 gives the line 
real power flows for a most heavily loaded single line ■ 
outage in IEEE 57 bus test system using the proposed methods 
as well as FBLP. The exact results obtained by NM load flow 
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are also given for comparison purposes. Table 3.5 shows the 
average error for the overall study and the highest of the 
maximum error observed for the different schemes. Besides, 
the relative CPU time in # is also given. In the analysis of 
average and maximum error, lines with less than 50 MW flow 
were not considered. The base case solution was used to 
start the outage calculation for the proposed as well as EDLP 
methods. Based on the representative results given in Tables 
3.4 and 3.5, the following observations are made. 

1. Overall simulation time of the proposed (I0 y , lE^) scheme 
is slightly faster than the (10, IV) scheme of ED1P. Simula- 
tion times are considerably reduced with the proposed (20 y , lE-^ 
scheme compared to the (20, 2V) scheme. The percent reduction 
in computer time becomes greater as the number of buses in the 
system increases. Although the average error obtained by us- 
ing the (20^, lE^) scheme is slightly higher than the (20, 2V) 
scheme, it is still less than the acceptable limit of 1 per- 
cent. This shows that fast and reasonably accurate solutions 
can be obtained by the proposed HEEDS, method. 

2. The proposed EDLR technique takes about 10 to 20 percent 
less computer time compared to ED UP. The average and maximum 
errors in line real power flows for the overall study using 
the proposed (If, 1e) scheme are comparable to or better than 
those obtained from (1®, IV) scheme of EDEP . These values are 
slightly more for the proposed (2f , 2e) scheme compared to the 
(20, 2V) scheme of EDLP. 
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3.6 CONCLUSION 

In this chapter, we have presented two methods for 
contingency evaluation. Method I is based on eliminating 
the load buses using sparse matrix techniques. Ibis method 
has definite savings in computation time compared to PDLP 
if the number of P-V buses in the system is less than about 
10 f<> of the total number of buses. It requires more memory 
storage than PDLP due to the requirement of the storing of 
factors needed for updating the matrix [B^J . However, this 
additional storage is not a limiting factor, since sparsity 
is exploited. Method II is based on the fast decoupled 
technique in. rectangular coordinates. In comparison with 
PDLP, it requires 10 to 20^ less computer time. The core 
storage requirement is the same as in PULP. Network outage 
cases are simulated using the proposed methods in one or two 
iterations with an average error in real power flows less 
than 1 i° % These methods have been tested on the IEEE 14 bus, 
30 bus and 57 bus systems on DEC 1090 system. 



CHAPTER IN 


OPTIMAL LOAD PLOW SOLUTIONS AND 
ECONOMIC DISPATCH 


4.1 INTRODUCTION 

The optimal load flow (OLE) is a constrained nonlinear 
programming problem. The OLE and the simplified economic 
dispatch problems have been discussed in Chapter I (Section 
1.3). The methods reported for OLE are based on Kuhn-Tucker 
theorem [64], Lagrange multipliers [65,71], Hessian matrix 
[681 and reduced Hessian matrix 1 74-1 possessing second 
order characteristics. The various methods [80-83] proposed 
for economic dispatch mainly differ in the calculation of 
penalty factors. A brief review of these methods pertaining 
to OLE and economic dispatch has been made respectively in 
Sections 1.5.3 and 1.5.4. 

The objective of this chapter is twofold. The first 
is to utilize the NRM load flow technique developed in 
Chapter II for OLE solution in the method of Dommel and 
Tinney [65). This approach avoids the re-tri angulation of 
the Jacobian matrix at each load flow iteration and hence 
considerable saving in computer time is achieved. Then, the 
NETR load flow technique developed in Chapter II is used for 
economic dispatch. The optimization is carried out on an 
equivalent system consisting of only P-7 buses. It uses the 
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constant factors of the Jacobian matrix relating to the 
reduced system. Therefore, it is fast, and competitive to 
the existing methods [80,82] using TDLP and is very attrac- 
tive for real time applications. 

The proposed methods have been tested on the adapted 
IEEE 14 bus, 30 bus and 57 bus systems and the results are 
compared with those obtained by Pommel and Tinney’ s method 
[651 in terms of solution accuracy and overall computing 
times. 

4.2 STATEMENT OE THE PROBLEM 


The variables and parameters of the power system have been 
separated into unspecified and specified quantities represen- 
ted by [x] and [y] respectively, in Chapter II (Section 2.2,1). 
Some of the specified variables [y] are adjustable. Hence the 

A ~[ 

vector [y ] can be partitioned as [y] = t i where [u] is the 

jjPlj 

control vector* and [pj is the vector of fixed parameters. 

Thus, the power system variables and parameters are 
grouped into three categories as summarised below : 


1. State vector or dependent vector [x] . This consists of 
the voltage angles of P-7 and P-Q buses and the voltage mag- 
nitude of P-Q buses. 

2. Control vector or independent vector [ul. This consti- 
tutes the adjustable real power generation available for op- 
timum dispatch, voltage magnitude of generator buses and 
transformer tap ratios. 
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3. Constant or specified vector [ p) . This consists of the 
constant load power, specified generator reactive power aid 
other constant network parameters. 

The general OLF problem can be stated as a mathematical 
programming problem given b y. 

Minimize the scalar objective function, 

3?(f xl , { u}) subject to 

i) Equality constraints imposed upon the system by load 
flow equations (vide Chapter II, Section 2.2.1) 

Jg([x] , t u], [ p]) J =0 (4.1) 

ii) Inequality constraints on some parameters of the form 



The objective function may be the overall cost of pro- 
duction or total system losses. For minimizing the cost of 
production, the objective function takes the form, 

n 

P 

P(txj, (Ul) = l (a o3 + b oj P S3 + o oj (Igj) 2 ) (4.4) 

3=1 
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For optimal reactive power flow, the objective function is 
the total system loss and is given by 

E([x3 , [ u]) = P Gs (Cxl , { u] ) (4.5) 

where Pq s refers to the slack bus real power generation. 

4.3 OPTIMAL LOAD PLOW SOLUTION - 3Y NEWT ON-RI OHAEDSON TECHNIQUE 

The OLE problem stated in the preceding section can be 
solved efficiently by Dommel and Tinney's method 1 651 . In 
this method, the load flow equations are solved by NM, which 
requires the calculation of the table of factors of the 
Jacobian at each load flow iteration. In this section, the 
application of NBM (vide Chapter II, Section 2.5) is proposed 
in the place of NM. The proposed method differs in two as- 
pects from the method of Dommel and Tinney. 

1) The table of factors of the Jacobian obtained at the 
end of initial load flow are treated as constant during the 
subsequent load flows. 

2) The equations pertaining to the load flow and Lagrange 
multipliers are solved by using NRM. 

4.3.1 Method of solution 

Firstly, the inequality constraints are not considered. 
Then, the Lagrangian function is formed by introducing as 
many auxiliary variables (otherwise known as Lagrange 
multipliers) as there are equality constraints in (4.1) and it 
is given as 
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L = F([x], [u]) +JX] 1 [g(tx], [u] , [p] ) ] (4.6) 

where [ X ] = vector of Lagrange multipliers. 

Optimality criterion : 


The necessary conditions for the minimum of 1?, while 
satisfying the equality constraints [ g] are. 


tl\l- 

1 ax J 

= [g([xl 

» C ul, 

CP! ) ! 

= 0 

(4.7) 

[ L X 1 = 

1 3 x J 

— r 

1 3x j 

+ [ 

1 9x J 

T 

CX ] 

= 0 

(4.8) 

[L u ] = 

l a u 1 

= t 9 ”l 
l a u J 

+ [ 3 -£ 
l a u 

T 

) C X] 

= 0 

(4.9) 

The set 

of non 

-linear 

simultaneous 

equations 

(4.7 to 4.9) Is 


solved by iterative solution technique to obtain the optimal 
solutions. 

4.3.2 ITEM solution procedure 

The solution algorithm proceeds as follows : 

■Step 1 : Start the iteration count for the gradient step 
. k = 1 

Step 2 : Read the system data. Assume a starting value for 
the control vector [u] and for the unspecified 
voltage phase angles and magnitudes. 

Step 3 : Check if k > 1. If it is, proceed to step (5); 

otherwise continue to the next step. 



88 


Step 4 : Solve the load flow equations (4.7) by ITEM for the 
state vector Cxi. Obtain the table of factors of 

y v 

the Jacobian matrix corresponding to the coverged 
solution and store it. Proceed to step (6). 

Step 5 : Solve the load flow equations in (4.7) by NHM 

making use of the table of factors available from 
step (4). Proceed to step (7). 

Step 6 : Solve the equation (4.8) for [X] by a repeat solution 
using the table of factors available from step (4). 
Proceed to step (8). 

Step 7 : Solve the equation (4.8) for [X] by ITEM making use 
of the factors available from step (4). This re- 
quires 2 repeat solutions for 37BM - scheme I, 

Step 8 : At this point ([xj, [uj , [\1 ), compute the reduced 
gradient vector [1^1 from equation (4.9). 

Step 9 : If [1^1 is sufficiently small, indicating that no 
further improvement in the objective function is 
possible through the control of [u } and therefore 
the minimum is reached. At the minimum, the gra- 
dient vector components will be 


1 . = 
Ul 

0 

if 

uf ] 1 * 

u i < 

Xi . x 
Ul x 

0 

if 

u i = 

u max 

X 

I . i 
Ul ' 

0 

if 

u i = 

u min. 


Otherwise, proceed to the next step 
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Step 10 : Adjust the control vector u In the negative 
gradient direction using the formula 

Iu k+1 ] = [u k ] + [ a u k ) (4.10) 

[ Au k ] = - & k [l k ] (4.11) 

Advance the iteration count k = k + 1 and return 
to step (3). 8 is a scalar and is known as step 

length. 

4.3.3 Determination of step length 8 

The critical part of the solution algorithm given above 
is step (10). The computational efficiency and convergence 
characteristics are sensitive to the determination of the step 
length in the negative gradient direction. The step lengths 
are calculated at each succeeding, grad lent steps k by using a 
steepest descent algorithm f 721 as described below. 

Ir IT— *1 

Step 1 : Compare the signs of and . If they are equal, 
then 8 k = 1 .5 8 k ~^ . 

Step 2 : Otherwise, compute the inverse Hessian matrix diagonal 

A 

element h^ to interpolate the optimum value which has 
been overshot due to the reversed direction of the 
gradient component 


. = (u k_1 - uF)/(L ! . 
11 v 1 i" v ui 


k-1 


- L k .) 
ui 
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(b) Functional inequality constraints 


These are inequality constraints except those on the 
independent or control variables and are taken into account 
by the penalty function method. In this method, the objective 
function is augmented by a function of the violated constraints, 

i.e. 


Fdx] , ful) = F (lx] , lu]) + \ R. 

j e m J 


(4.13) 


■where P q is the original objective function; m is the set of 
all violated constraints and R. is the penalty for each of 

ti 

the violated functional constraints 


E o - 

V h j - 

,maxx2 

h 3 5 

if 

h 3 > 

, max 
h d 

= 

Y. (h min 
3 J 

- v 2 

if 

h 3 < 

, min 
h 3 

where Y 

is the 

penalty 

factor. 



(4.14) 


The main difficulty in this method is the selection of 
suitable values for the penalty factors which influences the 
convergence characteristics and solution accuracy. To alle- 
viate this difficulty, a variable penalty factor approach as 
indicated in reference -t 741 can be used. In this approach 
Y. may be considered as a function of the iteration count. 

To start with, a small value is chosen and then successively 
larger values are employed and this gives more weight to 
the objective function in the first few iterations. The 
larger value of the penalty factor used, as the solution is 
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approached, increases the weight on the minimization of 
the penalty function. 

4.4 ECONOMIC LOAD DISPATCH SOLUTION BY NETWORK REDUCTION 
TECHNIQUE 

The load flow solution by network reduction technique 
(NETR) developed in Chapter II (Section 2.6) is used for 
fast economic load dispatch. 

4.4.'1 Economic load dispatch 

Here, the objective function is the production cost as 
given by the equation (4.4). The control vector [u] consists 
of only controllable real power generations. Other control 
variables are assumed to be fixed. 

Problem statement 

Determine [u ] (real power generation) which minimizes 

ir 

the objective function E( [x] , [u] ) subject to 

i) load flow equation 

tg([x] , f u], [pj )] =0 

ii) limits on control variables [u ] 

Ir 

The functional inequality constraints are not considered. 
4.4.2 Method of solution 

The real and reactive power balance equations (load flow 
equations) are solved for the initial assumed scheduling of 
real power generation and the converged solution is obtained. 
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If only, the real power scheduling is changed, the load bus 
voltage magnitudes will remain fixed, because of the follow- 
ing facts : 

i) The voltage magnitude variations at the load buses are 
not too sensitive to the -real power changes. 

ii) The reactive power equations depend weakly on the phase 
angles C © 1 of the voltages 

Hence, the load buses are eliminated by converting the 
loads into constant shunt admittances based on the converged 
voltage solution obtained, using the initial scheduling of 
real power generation and an equivalent system consisting of 
only P-V buses is obtained. 

The problem is now reduced to : 

Determine [ u ] which minimizes P r ([xJ, {u] ) subject to 

i) load flow equations 

[g r ([x],[ul, fe])J =0 (4.15) 

ii) limits on control variables [ u ) 

where the subscript r refers to the reduced system consisting 
of only P-V buses. 

Introducing the Lagrange multiplier vector C^ r J, the 
constrained minimization problem is transformed into the 
minimization problem with respect to control vector [ Up ] of 
the unconstrained Lagrangian function L . 



The necessary conditions for the minimum are 


' E rxr.l 

= £ 


[uj, lp| )1 

= £0] 

(4.17) 


= £ 

9L 9 E 

a - *-» Ji + 

9g r I 

^ *9 x 3 t x rl 

= £0] 

(4.18) 

'W 

= C 

9 L 9 F 

n' = ( 8 u ?] + 1 

P P 

9u p 3 {X r 3 

= £0 3 

(4.19) 


T he set of nonlinear equations (4.17 to 4.19) relating to 
an equivalent system consisting of only P-V buses ig solved 
iteratively to obtain the optimal solutions. 

4.4.3 KETR solution procedure 

The main steps of the solution algorithm are as follows : 

Step 1 : Start the iteration count for the gradient step k = 1 
Step 2 : Read the system data. Assume a starting value for 

the control vector £u J and for the unspecified vol- 

ir 

tage phase angles and magnitudes. 

Step 3 : Check if k > 1 ; if it is go to step (5); otherwise 
continue. 

Step 4 : Solve the load flow equation (4.7) by RETR load flow 

method, and obtain the converged solution. Store the 

matrices [7 ‘q.q.j (non-triangularized) and [ Y| T[ ] (tri- 
\ 3^7? 9p 

angularized) . Eorm [ *-■ ] = [ — J and store its 

° x 36 v 

table of factors. Go to step (6). 
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Step 5 : Solve the load flow equations by NRM using the 

3 iE? 

factors of { available from step (4). Go to 

step (7). 

Step 6 : Solve the equation (4.18) for [ X^jby a repeat 

solution using the factors of { j - ] available 
from step (4). Go to step (8). 

Step 7 : Solve the equation (4.18) for [X J by NSM using the 

T ' . ...... 

factors of t ’-^1 available from step (4). 

Step 8 : At this point ([x] , [u] , fX r j), compute the reduced 
gradient vector (I J rU pl from equation (4.19). 

Step . 9 : If [ is sufficiently small, minimum has been 

reached and go to step (11). Otherwise, continue to 
the next step. 

Step 10: Set the new values for u . as > _ ' 


r 

i 

max . „ 

u . ; if 

pi 

k 
u . 
Pi 

> 

max 

V 

k+1 _(v 
pi : 

"pf » . “ 

k 

V 

. k . 

+ A Upi < 

u m ^ n 

pi 

! 

( u pi + au pi 

) 

otherwise. 





, k fl k T k 

where a u pi = - .. 

Advance the iteration count k = k+1 and return to 
step (3). 

Step 11: Determine the load bus voltages from 

C Y£j.J f E L J = - fY LG I CEqI (vide equation 2.53). 
Now we have got the complete set of results. 
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4.5 NUMERICAL RESULTS OR TEST SYSTEMS AND COMP ARISON 

The adapted IEEE 14 bus, 30 bus and 57 bus test systems 
were used for obtaining OLE and economic dispatch solutions 
using the proposed NEM and UHER solution procedures respecti- 
vely and the results are compiled in Tables 4.1 to 4.4. The 
different test system specifications including the generation 
cost are given in Appendix D. The optimal solutions obtained, 
using the original Dommel and Tinney's approach based on MM 
load flow are also included. A comparison of the results 
listed in Tables 4.1 to 4.4 reveals : 

A. MEM Algorithm 

1) The total computing time is less by about 50 to 60$ 
as compared to by using MM. 

2) The total number of gradient steps, load flow iterations 
and the optimum points reached are practically the same as 
compared to those obtained by using MM. 

B. UETR Algorithm 

1) The overall computer time is about 10 to 20$ of that 
required by using MM. 

2) The saving in computer time as obtained in (1) increases 
as the system size increases and as the number of P-V buses 
in the system decreases. 
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TABLE 4.1 OPTIMAL LOAD PLOW RESULTS USING- NM, NRM 
AMD NETR ALGORITHMS 

IEEE 14 BUS SYSTEM 


\ Initial 
j schedule 


; _ , .. Op tima l schedule 

| NM \ NRM i 
1 algorithm i algorithm 5 


NETR al- 
gorithm 
(only real 
power 
schedule) 


ireal power 
\ generation 
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; p g2 
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\ 

l 
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1 69,44 
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1 
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1.06 

! V 2 

E 

1 

1.045 j 1.044 

1 

1.044 

; 1.045 
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\ 
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■? 
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I 
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* 
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= 

i 
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i 

> 
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i 
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i 
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X 

£ 
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; total gen era- 

5 

% 
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* 
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; 

i 
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1 
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I ARTE 4.2 OPTIMAL LOAD FLOW RESULTS USING- M, ITEM 11® 
ITETR ALGORITHMS 

IEEE 30 BUS SYSTEM 


(Initial 

jschedule 
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; algorithm 
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TABLE 4.3 OPTIMAL LOAD PLOW RESULTS USING MM, NRM 
ANL NETR ALGORITHMS 

IEEE 57 BUS SYSTEM 
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TABLE 4.4 COMPUT AT ION AL COMPARISON BETWEEN NM, NRM AND NETR 
ALGORITHMS FOR OPTIMAL LOAD FLOW 

(IEEE 14 BUS, 30 BUS AND 57 BUS TEST SYSTEMS) 
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3) The total number of gradient steps, load flow iterations, 
and the optimum values for real power generation are slightly 
different from those obtained by using MM, while the cost of 
generation of the optimum point is in close agreement. 

4) The assumption that the different real power schedules do 
not change the load bus voltages very much is quite valid. 

4.6 OOHCLUSIOW 

In this chapter, we have presented new solution pro- 
cedures for Old? and economic dispatch calculations based on 
MEM and METR load flow techniques respectively. It is 

t 

reported in reference [72] that the method using EDDP gives 
an overall advantage in speed of 2 to 1 over the method of 
Domrnel and Tinney. Since, the proposed MEM algorithm gives 

I 

an overall savings in computing times of more than 50# 
compared to the method of Domrnel and Tinney, it is competitive 
to the method using EDDP in terms of speed, while the Jacobian 
is available for further applications. The proposed METR 
algorithm for economic dispatch requires about 10 to 20# of 
the overall computer time compared to the method of Domrnel 
and Tinney. We may thus conclude from above, that this al- 
gorithm is atleast twice as fast as the one due to EDIT 172], 
and is thus ideally suited for real-time applications. The 
effectiveness of the proposed techniques is demonstrated on 
the adapted IEEE 14 bus, 30 bus and 57 bus test systems. 



CHAPTER Y 


TWO LEVEL SOLUTIONS EOR LARGE SCALE 
POWER SYSTEMS 

5.1 INTRODUCTION 

The large scale nature of the modern power systems 
requires the application of sparse -matrix arid/ or decom--- 
position techniques. A number of decomposition techniques 
are available in the literature which can be implemented 
in-, single [11, 27,33} or multicomputer [11,85 -87,96,971 con- 
figurations. In the last few years, there has been a consi- 
derable interest in hierarchical control methods in which 
the overall problem is split into a set of smaller sub- 
problems coordinated by a central agency. Brameller and Lo 
[ 92] proposed a diakoptical optimisation procedure in 
formulation using the method of Dommel and Tinney [65] . 
Although it is reported that the method is suitable for multi 
computer configuration, test results are not given in this 
reference. Sasson [90] presented a decomposition technique 
to solve the load flow problem using nonlinear programming- -• 
method. This method overcomes the convergence difficulties 
encountered in the application of the nonlinear programming 
load flow method to large scale systems. Happ proposed a d.c 
load flow [87] and real power dispatch [86] procedure in the 
multiarea case. This method of real power dispatch utilizes 
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the well known B-coefficients . JT.arita and Hammam [98] and 
Kobayashi et al [93] discussed the voltage and reactive power 
control for real-time applications amenable to multicomputer 
computation. More recently, Irving et al [107] made a com- 
parative study of the application of different linear pro- 
gramming methods for on-line real power dispatch using the 
hierarchical control theory. Bielli et al [108] proposed a 
two level technique for the economic real power dispatch 
based on the hierarchical approach. These methods by Happ, 
Irving and Bielli do not take care of the limits on the reac- 
tive power and bus voltages. 

The main objective of this chapter is to present a new 
two level optimal load flow procedure using the nonlinear 
programming method of Dommel and Tinney [65]. It optimizes 
both the real and reactive power using the adjustable real 
power generations,voltage magnitude of generators and trans- 
former tap ratios as control variables. Before describing 
the proposed optimal load flow method, the two level NM, ITEM 
and BDItP load flow solution algorithms are presented. The 
proposed two level technique combines the advantages of both 
decomposition and sparsity and is suitable for multicomputer 
conf iguration. The information relating to only the cut buses 
is transmitted between the two levels. The convergence cha- 
racteristics depend on the value of the acceleration factor 
used and the number of areas into which the original system 
is decomposed. 
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To implement the two level technique, a large scale power 
system is tom into several subsystems or areas by removing 
the connecting lines. Each area has its own dispatch computer, 
and in addition, there is a pool dispatch computer. There is 
no direct communication between the areas. The decoupling of 
different areas is achieved by introducing pseudo variables. 

In the first level the area solutions are obtained. The 
second level (pool dispatch computer) coordinates the areas :~ 
via ' pseudo variables. This process is iterated until the 
specified convergence criteria is satisfied. The proposed two 
level algorithms for the load flow and optimal load flow pro- v 
blems have been tested on the IEEE 30 bus and 57 bus systems, 
decomposed into 2 or 3 areas. 

5 . 2 LOAD PLOW 

5.2.1 Development of multiarea load flow equations 

The large scale power system is torn into N areas by 
removing the cut lines. The following assumptions are made for 
the system tearing. 

1 ) No mutual coupling should exist between the lines 
belonging to the different areas. 

2) Each area is individually connected. 

3) A common bus (usually the ground) is available for all 
the areas. 

The general guide for the system tearing is as follows : 
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1) Separation of the networks of different utilities, 
identifiable power pools etc. 

2) Employment of the various proposed network tearing 
algorithms £ 91 , 95 , 102 ] 


Eor the ith area, we define 

tv 1 ] 
te 1 ] 


IpJ] 

tQj] 

f o£] 

lx 1 ] 

lu 1 3 

i 

[P 3 

t p Li 


lw 0 1 


: bus voltage magnitude vector 
• bus voltage angle vector 
: real power generation vector 
: reactive power generation vector 
= real power load vector 
: reactive power load vector 

= state vector or dependent vector 5 ^Sectic^f^ 2) 
control vector or independent vector 3 
= parameter vector 5 

: number of generators 

= real power injection vector due to the cut lines 

= reactive power injection vector due to the cut lines 

n 3 s 


U— J 




= = cut line real power and reactive power flow 
vector 


(a) Gut line flow equations 


Phe real and reactive cut line power flows are treated as 
additional problem variables. They can be written -as a 
function of bus voltages and cut line parameters. The real and 
reactive power flows in the ith cut line interconnecting the 
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pth "bus in the area i and qth hus in the area 3 are 


c(p,Q) 


" V P ° 0S(9 P _ e 9 ) + b i Sifl(e p 




+ v 1 v 1 g + v 1 v 1 

■n -n o a ^ *r 


P P 2 


p p 


gl 

*2 

T~ 


(5.1) 


Q 


c(p,q) 


•where 


-Tp siP(»p - 0 3 q) 

-V 1 Y 1 b„ - Y 1 Y 1 ■*=■ 
p p A p P 2 ^ 


h A cos ( 9 ^ - Oq)) 


(5.2) 


y = + 3 bj = admittance of the 2. th cut line 
y> = g' + j h = total line charging admittance 

x* jL 

(b) Multiarea load flow equations 

The real and reactive power injections due to the cut 
lines in the area i are given by 


i p y - [0 p> ' p =i 

fe L> = tc ^'V 


(5.3) 

(5.4) 


where C C 1 ! = bus connection matrix which gives the connection 
of .out lines to the buses in area i* 


Ihe element C . (row a; column 3) is given as 
a J 

0 . = 1, if the cut line 3 is connected to the bus a 
a 3 

= 0, if the cut line 3 is not connected to the bus a 

The matrices [ 0 * ] and [ 0 ^ ] are the same only if there are no 
P-V buses present in the area i. Otherwise { ] is obtained 
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from [Cpl "by deleting the rows corresponding to the P-V 
"buses. 

Then, the real and reactive load flow equations for 


the area i will "be 

- [P^ehv 1 )] - {Cp P 0 i =o (5.5) 

-I« i (e i ,7 1 )] - [Cq][« 0 ] =0 (5.6) 

In concise form, the multiarea load flow equations can "be written as 
[g 1 ([x 1 ] , [u] 1 ], [p 1 ])] - [C^CWq] =0 (5.7) 

The cut line flow equations are given by 

[S c (l x1 ]> [ y1 ]» in, [u 2 ] ... (/], [if 1 ])] -[ V C 1 = 0 (5.8) 

Combining (5.7) and (5.8), we get 

t g 1 ([x 1 ] , {u 1 ] , [p 1 ] ) 3 - [C 1 ] tg c ([x ] , [u ] , {X 3 , [U ] ... 

[x 1 ], [u 1 ], [x W ] , tu F l )] = 0; i = 1,2, ... IT (5.9) 


As such, equation (5.9) cannot be solved independently , 
because of the coupling terms. We introduce psoudo- 
variables [<jl as per the following relationship to uncouple 
the terms in equation (5.9) 

[ah =[x i 3 (5.10) 

.A. 

loM = [ U 1 ]; i = 1 ,2, . .. IT > (5.11) 

Then, equation (5.9) reduces to 
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fg 1 (Cx 1 ), [u 1 }, [p 1 ] )] - ,[a5 , ... 

[x 1 !, [u 1 }, . ..[ oJ,[o5)l J i = ^2, ... IT (5.12) 

(c) First level 

The pseudo variables to^I , [<r^J (j = 1,2, ... IT, 3 ^ i) 
are assumed to he known as far as the area i is concerned. 

Then, the load flow equations ( 5 . 12 ) are solved independently 
for each area. 

(d) Second level 

The pseudo variables [o^] ,[ a*?] , (j = 1,2, ... IT) are 
determined using equations (5.10) and (5.11). 

(e) Acceleration Method 

To achieve better convergence for the iterative process, 
an acceleration method - Successive Over Relaxation (SOR) can 
be used in the first level or in the second level. However, 
the acceleration at the first level is found to decrease the 
number of iterations for convergence. 

Acceleration at the first level : 

At every calculation of [x 1 ] , a correction is applied in ? 
the following manner 


i k+1 i /, 1 k+1 

I A x 3 =0 ({ xl 


r k x 

- [xl ) 


(5.13) 


where a 1 is the acceleration factor for area i, with the value 
ranging between 1 and 2, and k is the iteration number. Then, 
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= [x 1 ] 151 +[Ax i ] k+1 ; i = 1,2, ... IT (5.14) 

Fig. (5.1) shows the transfer of information between the first 
and second levels for the load flow solution. 


5.2.2 Solution of multiarea load flow equations 

(a) Newton-Raphson Technique 

linearizing the equations given by (5.12) using Taylor 
series expansion and simplifying we get, 


i » 

; t ap 1 ] \ 


[ j i 


! [A e H \ 


(5.15) 


: [ aq 1 ] : 


i A Y i 

; [ ^ \ 

L. Y 1 


i=1 ,2,... IT 


where 


< 

S Ai 

If J 1 1 


= f 


r ji i i 

l J 3 =*- 


t^3 


3(P+P in }1 
38 


] 




* 

3(Q + a ) i I 

t- ) < 

* 


f J 2l 


= [ 


3(P + P in ) 1 j 

vr v i - 


38 


f JaI 


s(Q+Q ill ) 1 J 

f 3 V 


[AP 1 J = real power mismatch vector 

fAQ^-J _ reactive power mismatch vector 

[A0 1 ] = voltage phase angle increment vector 


{AY 1 ] = bus voltage magnitude increment vector 




3 1=1 to 4 


= submabrices of the Jacobian 
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(b) Eewt on -Richard son Technique 

Let t?<°n = Jacobian matrix evaluated using the 

initial estimates. 

Then applying Rewt on -Richard son technique (vide Chapter II, 
Section 2.5.1) for each subsystem, -we get the final converged 
value of [Ax j x ' at the kth major step as 


[Ax 1 ] ^ = [Ax 1 ] ^ + ... (5.16) 

i = 1,2, ... IT 

■where __ 

i i Tap 1 ! 

[«xM k) = ch%i > 

° [ a qj 

”** r 

r A *i] 

IJ^ « = ijt 051 ! - " - tj (1) ] (k) > 

' 1 |A C^j 0 

and so on. 

(c) East Decoupled Technique 

The coupling matrices I 1 and [ ] are neglected in 

equation (5.15) by application of the decoupling principle. 
The resulting equations are further modified by incorporating 
the suggestions proposed in the reference [29] to yield. 


t - 4 / 1 ] 


- t^ 1 ) Ue 1 ) 


(5.17) 


A Q 1 

~T~ J 


= fB” 1 } 


[AV 1 j 


= 1,2, 


IT 


(5.18) 
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where tB^ 1 ) = [B' 1 } + [Bi^l 
IB” 1 } =[B !,i ] + [B*£-J 

[ B l ^ } and [ B”^ } are the diagonal matrices with 
elements composed of only the cut line suscept ances. Since 
[B' 1 ! and [ B t,x I contain only the branch susceptances of the 
area i, the matrices [B^ 1 ] and [B^ 1 ] are real, sparse and 
constant and hence they need to he triangularised only once, 
at the start of the iterative process. 

5.3 OPTIMAL LOAD PLOW 

5.3.1 Multiarea power system optimization problem 

The multiarea optimal load flow problem can be stated 
as follows : 

Minimize the scalar objective function P of the system 
variables 

P = P 1 ([x 1 ] , [u 1 ! ) + P 2 ({x 2 }, [u2j ) + ... ) 

(5.19) 

subject to the following equality and inequality constraints, 
a) Equality constraints 
(1) Area load flow equations 

[g^x 1 ] , [u 1 ] , [P 1 ])] -[ wjj = 0; i = 1,2, ... N 


( 5 . 20 ), 
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(2) Area equivalent cut line injection equations 

[SqCEx 1 ] , [u 1 ] jtx 2 ] , [U 2 J ... [X F J, [U N ] )] -[ ¥^J = 0 

i = 1 , 2 , ... U ( 5 . 21 ) 

where 

l ¥j] = [ C 1 } [¥ c ] 

and 

[§c C lx 1 ! , t^ 1 ] , [x 2 ], [u 2 ]... [x F ],{U N ])] = 

[G 1 ] rg c (fx 1 ] , [ U 1 ] , [X 2 ] , [U 2 ] ... [X H J ,[U F ])j 

Ob) Inequality constraints 

These are the lower and upper limits on real and 
reactive power generation, voltage magnitude, transformer tap 
positions, slack bus power and line flows. In compact form, 
these are represented for the ith area as 


u^ . ^ 

3 min 


u. 




u 


3 max 


( 5 . 22 ) 


h 3 min 




3 max 


(5.23) 


The objective function is the overall cost of production and 

it takes the form for the ith area as 

n 1 


P 

3=1 


Phrx 3 ! , [U 1 ]) = V Cab + b * 3 p | 3 + 0 ?L) 2 ) (5.24) 


Using the classical optimization method of Lagrange - . multipliers 
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[651, the minimum of objective function F subject to 
the equality constraints is found by minimizing the aug- 
mented objective function known as the Lagrangian function. 

5.5.2 The Lagrangian function 


The Lagrangian function is formed as 
N F 



i=1 


[x F ], [ u iT ])] - [tfj] } (5.25) 

where 

[X 1 ] = vector of Lagrange multipliers for the ith area load 

flow equations. 

jXq ] = vector of Lagrange multipliers for the area equi- 
valent cut line injection equations. 

5.3.5 Two level formulation 

The Lagrangian function is decomposed by introducing the 
additional variables called as pseudo variables [0] as given 
by (5.10) and (5.11). This allows each area to be decoupled 
and hence solved independently . These equations result in the 
additional equality constraints represented by 

r a X ] - C X 1 J = 


0 


(5.26) 



and 
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K' - 


[u 1 ] 


= 0 


- 1 2 
_ ■ , , 


N 


(5.27) 


With the above equality constraints, the Lagrangian function 
1 Q is modified as 


L = 


L o + 


N 

I 

i=1 


it 


i xii 


T 


x 




tx 1 !}* l U i :.] 1 { [<$' - (U 1 ! } 

°-u 

i=1 


(5.28) 


where f X 1 ] and [X 1 ] are the Lagrange multiplier vectors 
a x 

for the equality constraints introduced by the pseudo varia- 
bles. 


The total Lagrangian is readily split into IT independent 
sub-Lagrangians by considering [ cr ^ ] andfcr^j as constant 
except for the area i as indicated below. 


IT 

L = l L i (5.29) 

i=1 

where 


I-i = ^(fe 1 !, Cu 1 !) + [X 1 ] 1 { [ g 1 ( [x 1 ] , [u 1 ], [P 1 ])] -tW^] } 


+ fx e 3 { fgc ([0 ^ »f a u J » f a xJ ’ f a uJ fx ^ » t ul 5 

( a x J > la § )l " [¥ S } ’ fu 1 I 

T-3* 

IT T IT 

+ l 1 X aJ p x3 + If X 1 Ha ^ J 
3=1 X 3=1 a u U 


D^i 


3^i 
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5.3.4 Optimality criterian 
(a) First level 

There are H areas in the first level of the hierarchical 
structure -with independent lagrangians. Necessary conditions 
for stationarity for 1 -j , L 2 , . .. are 

= - [^| + t^uW^t^-cx/^ 

x . 3 X 1 8 X 1 ax ax 1 c °X 

(5.30) 


i 3 L i 3F 1 T i 9g o i l 

t \t = t ~b = t ii - t ® i i» i+ 1- ti t\,i - 1\, i = 


(5.31 ) 


Hjl- t '1 =(g i ([x 1 l , tu 3 ) , [P 1 ) )l - 1 = 0 (5.32) 

ax 


i 3 1 . 

>¥ - t r i -(* l -i^J = o 

C 3 ¥ rt 

c 


(5.33) 


[L X J = E ax i' 1 =tgc ( f xij ^ C^ 1 ] »[4^f a u l)} - = 0 (5 ' 34 > 

0 o 

for i = 1 , 2 , . .., N; 3 = 1 , 2 , ... N and 3 £ i 
Substituting for[X^] and rlj obtained from (5.33) and 
(5.34) in equations (5.30 to 5 . 32 ), we get 


=i 


aCg^gi) * 4 i 

1 + c , ] r x 1 ]- cx i = 0 


(5.35) 


rT i a F 

[1 J = r i ] + t 

3 u 


a(g i + gi) T i i 

■ i 3 fx 1 - [X 1 = 0 

3u x a 


(5.36) 
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CIu I = Cg 1 (f x 1 ! , Cu 1 }, [p 1 ] )] - Cg^jx 1 ] ^[u 1 ] ,[a^] = 


0 


(5,37) 


for i = 1,2, ... H; 3 = 1,2, ... F and j / i 

For any given values of [ X 1 ] ,[X 1 ], and f 0 ^ j determined 

at the 2nd level, equations (5.35 to 5.37) give the corres- 
ponding optimal decisions [x 1 ] , [u 1 ],^ 1 ) and hence [ X^ ] at 
the first level. 


(h) Second level 


The total Lagrangian system equation (5.29) is optimized 
at the second level. Fecessapy conditions for stationarity 
for 1 are 


[ 1 

3X a 1 
a x 

= !^]- [x 1 ] 

= 0 


(5.38) 

31 . 
3X 1 1 
°u 

= t oj] - [ u 1 ] 

= 0 


(5.39) 

3a 1 J 

X 

F 

“O + J l 

a x 3=1 

dA 

3g c ,T 

8 a 1 1 

X 

t*cJ - 0 

(5.40) 

1 B<^ ] 

u 

for i = 

F 

=[X i ] + I [ 

O 

u 3==1 

3»4i 

1,2, ... F. 

3 «0 I 

I l o 1 = 0 

(5.41) 
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Por any given values of [X 1 ] , 
the first level, equations (5 
ponding optimal decisions [oa 


[u 1 ] and [xb determined at 
c 

.38 to 5.41) yield the corres- 

i . foil > t^„ 1 1 3114 tx„ x ) 


at the second level. 


5.3.5 Solution technique 

The results summarized by equations (5.35 to 5.41) are 
used to solve the two level optimization iteratively in three 
major steps as indicated below. 

Step (1) : Calculation of state vector {x J 

•Assume a starting value for the control vector [ u], 

Then, the state vector [x] is obtained by solving the multi- 
area load flow equations represented by (5.37) in the first 
level and the equations (5.38) and (5.39) in the second level 
iteratively in exactly the same manner as described in 
Section 5.2. 


Step (2) : Calculation of Lagrange multiplier vector [X] 
Prom (5.35) and (5.40) we get 


[ * 1 = - [ 


3(g X + gt) " 1 


c y T 3 P 1 . 1 

” r~ ] { [ . 1 I - i x a ] } 

3x 3 x CT x 


E* 


1 r 

a 1 - - .1. 

°x n=1 


3g c T j 

ti l ) E * 1 

3=‘ 

jjii x 


(5.42) 
. .. Pirst level 

(5.43) 


for i = 1,2, ... N 


. . Second level 



119 


For the converged value of [ x ] available from step (1), 
equations (5.42) and (5.43) are solved iteratively to obtain 
[XI. To speed up the convergence rate, an acceleration 
method (SOR) is used in the first level. 

■Step (3) : Calculation of control vector [ u] 

For the set of values [x ] and [ X ] available from 

steps (1) and (2), ^ ^ } is calculated at the second level 

tfu 

using equation (5.41). Then, in the first level [1^ ] is 
calculated using equation (5.36). If is sufficiently 

small, indicating that no further improvement in the objec- 
tive function is possible through the control of [u 2 ] and 
therefore the minimum is reached. Otherwise, adjust the 
control vector in the negative gradient direction using the 
formula 

=( u 1 ]* + [A u 1 ]* (5.44) 

■where 

[A u4 k (5.45) 

i = 1 ,2, . . . N 

and return to step (1). g is a scalar and is known as the 
step length. Its value can be chosen in a similar manner as 
indicated in the preceding chapter, Section 4.3.3. 

Fig. (5.2a-c) depicts the transfer of information between th 
first and second levels at the three major steps of the opti 
mization. 
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5.3.6 Inequality constraints 


The inequality constraints encountered in the optimiza- 
tion process are given as : 

a) Parameter inequality constraints 

These are the constraints applied to individual control, -• 
variables namely 


^ u min 1 ^ 


[- 1 ] 


*<ax 


1 


= 12 


N 


As indicated in Chapter IV, Section 4.3.4, in the event that 
the correction A u^ found from eqn. (5.45) causes the new 

t) 

value of the jth control variable in the area i, i.e. u^^ + ^ f 

d 

to exceed one of its limits, the solution algorithm sets u i 

d 

to the corresponding limit. 


u 


i(k+1 ) 


U 1 

. { 3 max’ 

= 4*- 

I 3 min’ 


if(u^ k ^ + 

d 

if( U i (k) + 

d 


aui (k) ) 

d 

Au^ k ^ ) 

d 


i u 


i(k) 




otherwise 


> uf 

3 max 

< u* 

3 min 


b) Functional inequality constraints 


The inequality constraints not involving the control 
variables are taken into account by the penalty function method 
as indicated in Chapter IV, Section 4.3.4. The objective fun- 
ction is augmented for the ith area as 



124 


where F q is the original objective function and m is the 
set of all the violated constraints in the area i and R^[ 
is the penalty function for each of the violated functional 
constraints in the area i 


Rt = 


*l< h J 


h 1 V 
j max ; 


if h: 


h 3 : 

3 max 


Y 1 (h 3 : . - h^)‘ 

3 3 mm 2 


if h. 


< hi 


3 mm 


where y 1 is the penalty factor and h 1 is the compact repre 
sentation for the functional inequality constraints in the 
area i. 

5.4 NUMERICAL RESULTS ON TEST SYSTEMS AND COMPARISON 


The proposed two level computational method was applied 
to the adapted IEEE 30 bus and 57 bus test systems shown in 
Figures 5.3 and 5.4 respectively. The system specifications 
are given in Appendix L. These systems were decomposed into 
two as well as three areas. The load flow results obtained 
by solving the multi-area load flow equations using NM, ERM 
and FELP techniques are shown in Tables 5.1 and 5.2. Tables 
5.1 and 5.2 show the total number of iterations, the relative 
CPU time in percentage and the optimum acceleration factor 
for the different systems studied. Table 5.3 gives the 
effect of acceleration factor on the number of load flow 
iterations needed for the convergence of IEEE 57 bus 
system. For comparison, the results obtained based on the 
computation without decomposition are also given. 




DECOMPOSED INTO 


3 AREAS 
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FIG- 5 -A IEEE 57-BUS SYSTEM DECOMPOSED INTO S 

3 AREAS. 1 



TABLE 5.1 COMP UT AT ION AL COMPARISON OP LOAD PLOW SOLUTION BETWEEN 

COMPUTATION WITHOUT DECOMPOSITION AND TWO LEVEL COMPUTATION 
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TABLE 5.3 EEEECT OE ACCELERATION PACT OR ON THE NUMBER 
OE ITERATIONS EOR CONVERGENCE 


IEEE 57 BUS SYSTEM SPLIT IN TWO AREAS 
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TA3LE 5.4 OPT MAE LOAE LOW RESULTS USING TWO LEVEL 
COMPUTATION 


IEEE 30 BUS SYSTEM 

System split in > System split in 

two areas i three areas 
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TABLE 5.5 OPTIMAL LOAD PLOW RESULTS USING TWO LEVEL 
COMPUTATION 

IEEE 57 BUS SYSTEM 
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TABLE 5.6 COMPUTATIONAL COM ARISON OP OPTIMAL LOAD PLOW SOLUTION BETWEEN 
COMPUTATION WITHOUT DECOMPOSITION AND TWO LEVEL COMPUTATION 

IEEE 30 BUS AND 57 BUS SYSTEMS 
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r o.'he optimal load flow results are given in fables 5*4 
to 5.6. The optimum values of the control variables together 
with generation cost are given in Tables 5.4 and 5*5. Table 
5.6 gives the total number of gradient steps and the relative 
CPU time in m is used for obtaining the load flow solution 
at each gradient step. In all the above calculations, acce- 
leration is used at the first level. The optimal load flow 
results obtained using computation without decomposition have 
been given in the preceding chapter (vide Tables 4.2 and 4.3). 

A study of the results presented in Tables 5.1 to 5*6 
reveals the following points. 

A. load flow 

1) The computational time per load flow iteration of the 
proposed method is faster than by the method without decom- 
position. 

2) The proposed method takes more number of iterations for 
the load flow solution compared to the method without decom- 
position, 

3) The computing time per iteration, total number of itera- 
tions and hence the overall computer time vary with the number 
of areas into which the system is decomposed. The convergence 
rate is adjustable and depends on the value of the accelera- 
tion factor used. The optimum acceleration factor changes 
slightly with the number of areas and the type of the load flow 
t echni que used. However, its value is found to be between 1.4 
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to 1.6 for the different systems studied. This value is 
similar to that used in Gauss-Seidel load flow. 

4) In general, the total computer time using HM, HEM and 
IDIiP load flow techniques is comparable to that of compu- 
tation without decomposition. 

B. Optimal load flow 

1) The results obtained agree very well compared to those 
using computation without decomposition. 

2) The number of gradient steps needed to reach the optimal 
solution is nearly the same as compared to the method without 
decomposition. 

3) The total number of X -iterations depends on the acce- 
leration factor used and the number of areas. Its optimum 
value is found to be between 1.4 to 1.6 as in the case of 
Gauss-Seidel load-flow. 

4) The overall computer time is comparable, when the "system 
is solved as a whole or decomposed. 

5.5 COFCLUSIOH 

A two level formulation is used for the analysis of load 
flow and optimal load flow in large scale power systems. The 
attractive features of the technique are that 

i) the optimization of both real and reactive power 
is possible. 
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ii ) the computational advantages of sparse matrix 

techniques are exploited in the analysis of areas, 

lii) there is a significant reduction in the computer storage 

since now, we have to store only the network matrices for 
each area 

iv) the data transmitted between the, 

eeween the area computers and pool 

computer is reduced. 

v) the overall computer time is comparable to that of 
computation without decomposition. 



CHAPTER 71 


CONVERGENCE CHARACTERISTICS OE THE NEW LOAD PLOW 

SOLUTION METHODS 


6.1 INTRODUCTION 

We have presented three load flow solution methods namely, 
EDLR, NRM and NETR in Chapter II. In this chapter, we study 
the convergence characteristics of the proposed methods to 
examine the existence and uniqueness of solutions. Hollowing 
the lines of Wu {33], the proposed load flow solution methods 
are represented individually in the form of a nonlinear itera- 
tive scheme. Then, the convergence conditions are derived 
based on the convergence theorem which has been obtained from 
the fixed point and generalized mean value theorems. Some ob- 
servations about these conditions of convergence and uniqueness 
of solutions for the proposed load flow methods and also EDLP 
are made based on the test results obtained on the IEEE 14 bus 
test system with different R/X ratio of a. transmission- line. 
Before deriving and verifying the convergence conditions for 
the proposed methods, a few existing studies available on the 
convergence analysis are briefly reviewed. 

6.2 REVIEW OE STUDIES ON CONVERGENCE ANALYSIS PROBLEM 

Meisel and Barnard [12] investigated the convergence 
characteristics of NM and modified NM based on Kantorovich 
theorem. Korsak 1 1 3j gave the definition and stability of the 
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load flow solutions and analysed the possibility of more than 
one stable solution. Wu [33] presented a theoretical study 
of the convergence characteristics of EDItP. He showed that 
the convergence conditions are dependent on R/X ratios of the 
transmission lines, scheduled power injections etc. Two al- 
ternate upper bounds for the contraction mapping constant, 
which is merely a number calculated to examine the convergence 
conditions have been derived. However, these bounds provide 
highly conservative estimate for the convergence region. Re- 
cently, Johnson [35] discussed the multiple and false solutions 
of the load flow problem. He reported that false solutions 
may occur due to a system having heavily loaded high voltage 
lines or generators feeding capacitive loads. Abe et al [ 38 3 
investigated theoretically the convergence region of NM in 
the polar and rectangular versions. The necessary conditions 
for HM to converge to a stable solution under certain opera- 
ting conditions are given. Successful convergence is reported 
for the cases having large phase angle corrections but not 
corrected in the 1st iteration. They have concluded- that 
the load flow convergence will be degraded if a starting value 
for the voltage less than the flat value is chosen. 

6.3 MLIUEAR ITERATIVE SCHEME 

As given in reference [33], we wish to solve a system of 
n equations f^ (x-j,x 2 , ... x ) = 2^; i -1,2, ... n in n un- 
knowns. In vector form, the system of equations are represented 
by 
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P(x) = Z (6.1) 

The suitable recursion formula to obtain the solution of 

A 

F(x) = Z is 

x k+1 =*(x k ) (6.2) 

1 2 

This generates a sequence of points (x ,x , ...) which con- 
verge to the solution point x s of the given system of equa- 
tions. Hence, at the solution point 

x s = $(x s ) (6.3) 

the solution of F(x) - Z = 0 must satisfy x - <{> (x) = 0. 

Hence, we consider the following relation, 

F(x) - Z = [ A(x)3 (x - <J> (x)) (6.4) 

where [ A(x)1 is a nonsingular matrix of order nxn and the 
elements of C A(x) 1 are functions of x. 

Equation (6.4) can be written as 

4 >(x) = x - [A(x)] _1 (F(x) - Z) (6.5) 

Thus, the general form of iterative scheme is given by 

(\r) 

x (k+ 1) = x (k) - tA(x) (k) 3 ” 1 CF(x) - Z) (6.6) 

Depending on the particular choice of the matrix { A(x)l , 
different algorithms are obtained. Generally, the load flow 
algorithms can be regarded as A nonlinear iterative scheme of 
the form given by (6.6). 

We express the proposed new 1 load flow algorithms into the 
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above general nonlinear iterative scheme as indicated in the 
following section. 

6.3.1 PDLR - as a nonlinear iterative method 


The final load flov; equations as given in Chapter II, 
Section 2.4 are 


l AP| 

= £ B 1 3 £&f 1 

(6.7) 

£e ] 

= (£V 2 1 - [f 2 J )* 

(6.8) 

£ aq J 

= £ B" ] C AeJ 

(6.9) 


The iterative algorithm is given by 
r 

If"! : [ B 1 ] 


; -pk+1 s 5 | ; r , i —1 4 

J i_„f 


| k+1 1 i k * 

i 0 J f 6 


1 r 

oli 


AP(f k ,e k ) 


-It 


Af^/4?k+1 k+1 kxr 

- AQ(f ,e v ,e )t 


1 t e j , 0 t [B” J 

X- ' Urn \ * 

Comparing equations (6.6) and (6.10), we get 


( 6 . 10 ) 


l ib'j ; o | 

t A(x) k j = j r j 

! o \ t b w j i 

L_ > 

In the evaluation of £AQ], newly evaluated [f J from equa- 
tion (6.7) and [e y ] from equation (6.8) are used. The matrix 
[A(x)j is constant and contain only the network susceptances. 

6.3.2 NRM - as a nonlinear iterative method 


Using equation (2.46) given in Chapter II, Section 2.5, the 
iterative algorithm is expressed as 
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f k+1 ! 

i 8 


fe k ! 


*1 - tJ^°h 


-1 


f“ ^ 

I - AP(e k , v k ) i 


\ yk+1 f 




l - AQ(e k , Y k ) i 


- [J (0) ] 


| j - AP(e k ,v k ) | 
-1<* h- : 


- AQ(8 k ,V k )| 


■ IJ (k) ] [J (0) ] 


-1 


) -AP(e k ,ir k )p 

; -;b 

1 

j u 

j -AQ(8 k ,V k )-5 




(6.11) 


X ♦ 0 # 


; k+1 J 

’ 9 ! 


bH 


a k i 
0 ; 


T k ! 




I -II 1 (9 k ,V k ) j 


AQ' (0 k ,Y k )_: 


Comparing equation (6.11) with (6.6), we obtain 


U(x) k ] =[J (0) ] 


The algorithm therefore, is a fixed Jacobian Newton algorithm 
with modified power injection equations. 

6.3.3 NETR - as a nonlinear iterative method 

The final load flow equations as given in Chapter II, 
Section 2.6 are 


AP v 

I ~y J 


° I B r 3 IA0 v 1 


(6.12) 




- 'hs» tE 0 > 


( 6 . 1 ?) 
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Equation (6.12) refers to an equivalent system consisting of 
only P-V buses. 

Equations (6.12) and (6.13) can be -written as 


f AP tt ] =[ V w . By CA0 Tr J 


(6.14) 


-1 


I 2 !! ~ f Y LG- H s q f 2 ! 1 (6.15) 

Using equations (6.14) and (6.15), the iterative algorithm 

can be written as 

i -r- 

k 


Se * +1 ! | 8 v| jtT, . Bp' 1 ! 0 !;.«(»*, y„) 


V ? 


V 


: ‘ 5 

1 l t 111 Jj \ \ 

—i L L 


i 


? r 

i i 




(6.16) 


Comparing (6.16) with (6.6) we get 

h.-fi »i 

[A(x) k ] = 5 1 "*" " r 

l 0 ! t*kiJ 


The iterative algorithm for the reduced system consisting of 
only P-V buses is given by 

te* +1 ] = te*3 - IT V . B^r 1 I- 4 P(B^, V v )I (6.17) 

Here, 

[ A(x) k l = I V v . B^l 
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6.4 CONVERGENCE CONDITIONS 

Using the convergence theorem defined in Appendix C, the 
conditions of convergence for the proposed load flow methods, 
are presented in this section. 

6.4.1 Convergence criterion for EDLR 

Suppose the region of interest for finding the solution 
to load flow problem with a flat voltage start (e^+jjf^ = 

1.0 + 30.0) is defined as 

|f ± | < e (6.18) 

| e^-1 |<e ; for i = 1,2, ... n (6.19) 

and C is the contraction constant given by 
c 

C c = max j |D <}> (f,e) [ j 

where 

; i j 0 ! • t b ']" 1 ; 0 j 

D <j>(f,e) =; I I - r t - f J(f , e) 3 (6.20) 

j 0 ! I l j 0 J IB"]” 1 ? 

J§* : 

Suppose that 
1) C 0 < 1 

ii) the initial correction Af?, A e? satisfying 
max { |Af? | , |AeJ |} < (1-C c ) e 


Then 
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1) there exists a unique solution within the region defined 
by (6.18 - 6 . 19 ) 

2) the iterates (6.7 to 6.9) converge to the solution. 

6.4.2 Convergence criterion for NBM 

Suppose the region of interest for finding the solution 
to load flow problem with a flat voltage start is defined by 

J ej < a for i = 1,2, ... n (6.21) 

I Y^-1 j < e for i = nw-1 ... n (6.22) 

and C is the contraction constant given by 

C c = max I |D 4 (0 ,V) f| 

where i- . — . 

i A t I 
c I i 0 •' 

• 1 ' > -1 

D *(e,v) = - cj ( 0) J [J(e,v)] (6.23) 

I } * ’ 

I 0 ! 1 ; 

i_ t j 

Suppose that , 

i) 0 o < 1 

ii) the initial correction AY? satisfying 

max { e/S} A9°| , | AY? |} < (1-C c ) e 

Then, 

1) there exists a unique solution 

within the region defined by (6.21 - 6.22). 
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2) the iterates (6.11) converge to the solution. 
6.4.3 Convergence criterion for HEIR 


Here, the conditions are derived for an equivalent system 
consisting of only P-Y buses. 

Suppose the region of interest for finding the solution 
to load flow problem with a flat voltage start is defined by 

jej < a for i = 1,2, ... n (6.24) 

I V ± - l| < e for i = m+1 ... n (6.25) 

and C is the contraction constant given by 

cr 

C cr = max || D 4 (0 y , \) | | 

where 

D4 (0 v ,V v ) = CXI -[ V v . B^ _1 tJ(8 v ,Y v ) 1 (6.26) 

and the subscript r refers to the reduced system. 

In the above equation, (Y^! the voltage magnitude vector 
relating to the P-Y buses is specified. 


i) 

ii) 


Suppose that 
C < 1 


cr 

the initial correction A9 


max 

i 


{ 


A. - 0 

e | A9vi| 


o 

vi 


> 


(1 - C 


a 


cr 


satisfying 

) t 


then, 
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1) there exists a unique solution within the region defined 
by (6.24 - 6.25) 

2) the iterates (6.12 and 6.13) converge to the solution. 

6.5 SYSTEM STUDIES 

The convergence conditions derived in the preceding 
section pertaining to the proposed load flow methods have been 
tested on the IEEE 14 "bus test system with different R/X 
ratio of ac- transmission line'. The maximum values of the 
contraction mapping constant C at any particular iteration 

C 

obtained using the proposed methods as well as EDEP are ta- 
bulated in Table 6.1. 

Discussion of the results 

The values of the contraction mapping constant 0 are 

c* 

less than 1 for all the proposed methods on the IEEE 14 bus 
test system when the line resistances are neglected. These 
values are more than 1 for ZDDR and less than or comparable to 1 
fan ; -TRM arid RET! when the line resistances are not neglected 
and R/X ratio of a heavily loaded line is increased on the 
IEEE 14 bus system. The 0 values obtained by using EDLR 
for the different cases studied are comparable to those 
obtained by using EDLP. Therefore, it is concluded that the 
convergence of EDLR is dependent on R/X ratios of the trans- 
mission lines as in EDIiP, whereas it is independent for SUM 
and NETR as in NM. That means, when R/X ratio is small, 
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TABLE 6.1 CONTRACTION CONSTANT FOR THE LOAD FLOW SOLUTION 
METHODS 

IEEE 14 BUS SYSTEM 


IMethod 

S 


i Contraction constant C - 1 

; o~ i 

| IEEE 14 Bus | IEEE 14 Bus j IEEE 14 Bus j IEEE 14 bus j 

\ system with \ system j system with j system with ; 

j all the re- ! \ R/Z ratio of I R/X ratio of I 

\ sistance 1 \ the heavily 1 the heavily \ 

j values ne- ] \ loaded line '\\ loaded line \ 

! glected j { increased by \ increased by? 

i « j 2 times i 5 times j 



•FDLP » 

J 

0.562 

’ 

1.945 } 

1.993 

2.354 

* i 

t j 

Ifdlr | 

3 

J j 

• 

0.744 

* 

j 

2.245 j 

j 

2.282 

2.764 

5 J 

!ra : 

\ * 

| 

0.585 ; 

i 

] 

: 0.793 | 

j 

0.805 

; 1.027 

i 

< 

* 

\ 

INETR i 

0.170 1 

' * 

» * 

0.209 | 

0.214 

: 0.397 
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faster convergence may be obtained by using FDLR, similar 
to FLLP, This fact 'has already been proved by the test re- 
sults. given in Chapter II, Section 2.7. For the systems 
having transmission lines with fairly large R/Z ratios, the 
larger values of C c for both FLIP and FDLR indicate that 
there may be convergence difficulties. This conclusion is 
similar to that reported for FDEP in reference C33]. 

6.6 CONCLUSION 

In this chapter, we have derived convergence conditions 
for the proposed FDLR, NRM and NETR load flow methods following 
Wu's approach [33] . The convergence of FDLR is degraded when 
R/Z ratios of the transmission lines are increased, which is 
similar to that for FLIP. The convergence of NRM or NETR is 
unaffected with higher R/Z ratios of the transmission lines 


as in NM 



CHAPTER VII 


CONCLUSION 


7.1 GENERAL 

This thesis has been aimed at developing improved solution 
methods for power system studies involving the load flow, con- 
tingency evaluation and optimal power dispatch and two level 
techniques for the load flow and optimal load flow based on 
decomposition and multicomputer conf igurations. 

We shall first review the significant results obtained 
due to the various suggested methods and then give the scope 
for further research. 

7.2 REVIEW OE THE SIGNIFICANT RESULTS 
Load flow : 

In Chapter II, three new load flow solution methods 
namely, EDLR, NRM and NETR have been proposed. 

i) EDLR has been derived using linear transformation from 
EDLP as well as direct from NM. 

ii) NRM has been developed using the simplified NM and 
Newton-iterative methods. 

iii) NETR has been developed based on network reduction 
techniques and NRM. 

Erom the test results, it is concluded that the above 
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methods have certain superior features compared to the 
existing widely used methods. These features are summarized 
below : 

i) EDIR is best suited for moderately accurate solutions 
compared to EDEP with 5 to 10$ overall saving in computer 
time. The memory storage requirements are similar to 
EDLP. 

ii) NEM is about two times faster than KM and with the 
availability of Jacobian, it is an attractive alternative 
where KM is preferred. 

iii) KEER is a comparable alternative to EDLP in terms of 
computational speed and core storage for medium size 
systems with convergence characteristics similar to KM. 

Contingency Evaluation : 

The proposed EDLR and KETR load flow techniques have been 
extended for contingency evaluation in Chapter III. In the 
EDLR algorithm, the matrix inversion lemma has been used 
whereas in the KETR algorithm, the matrix inversion lemma 
and matrix updating algorithm have been used to make use of the 
base case factors for the simulation of the different outage 
cases as in prevalent competing methods [29,57]. Digital 
simulation and comparison have shown that : 

i) EDLR is superior in computational speed by about 10 to 
20$ over EDLP. 
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ii) N^/IR is comparable or faster than PDLP in terms of speed, 
depending on the number of P-7 buses in the system. 

iii) the average error in the line real power flows is com- 
parable to EDIP for the different proposed iterative 
schemes using PDLR and NETR. 

Optimal load flow and economic dispatch : 

In Chapter 17, the 01P and economic dispatch problems 
have been solved. In the OLP problem, the equations pertaining 
to the load flow and lagrangian multipliers have been solved 
by using NRM. In the economic dispatch algorithm using NEIR, 
the optimization has been carried out on the reduced system 
consisting of only P-7 buses, based on the assumption that 
the real power schedules do not change significantly the load 
bus voltages. Numerical computation and comparison have shown 
that : 

i) there is an overall saving in computer time of more than 
50$ when NRM is used in the place of NM using the method 
of Dommel and Tinney {65] , without affecting the conver- 
gence characteristics, accuracy and general formulation 
of the problem. .Although, it is comparable to the method 
using ZDIiP [72 j in terms of computing speed, the Jacobian 
is available for further applications. 

ii) the approach by NETR for economic dispatch results in a 
great saving in computer time of about 80$ over the 
method of approach by NM [65] and hence is competitive to 
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the method using ILLP {73 . This shows the feasibility 
of this method for real-time applications. 

Two level solutions : 

In Chapter V, a two level optimization procedure has 
been developed for large scale power system optimization 
problem. In the proposed technique, both real and reactive 
power optimization have been carried out, considering the 
adjustable real power generations, voltage magnitude of 
generators and transformer tap ratios as control variables. 

This has been achieved by introducing pseudo variables. The 
original system has been split into a number of areas by the 
cut -lines. The area optimal solutions have been obtained at 
the first level and the coordinating pseudo variables and the 
associated Lagrange multipliers have been solved at the 
second level. Two level load flow procedures using KM, NHM 
and ILLP techniques have also been indicated. It is esta- 
blished that the proposed two level techniques give good 
results, for the load flow as well as OLF problems. The nu- 
merical results have shown that : 

i) the convergence characteristics of the two level technique 
are dependent on the acceleration factor used and the 
number of areas into which the original system is decom- 
posed. 

ii) in comparison with the computation without decomposition, 
the results obtained and the overall computer time are 
comparable. 
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Convergence analysis : 

A theoretical analysis of the convergence characteristics 
of the proposed FDLR, MM and 2IETR load flow iterative schemes, 
following the procedure indicated by ¥u (33], has been made in 
Chapter ¥1, to examine the existence and uniqueness of solu- 
tions in a specified region. The convergence conditions have 
been derived and verified on the IEEE 14 bus test system with 
different R/Z ratio of the heavily loaded line. The test 
results have shown that : 

i) there may be convergence difficulties encountered, with 
large R/Z ratios of the transmission lines using FDLR. 

This is similar to as reported for FDEP. 

ii) the convergence of FRM or EETR is not affected by the 
R/Z ratios of the transmission lines as in M. 

7.3 SCOPE EOR FURTHER RESEARCH 

In reference [110], a new algorithm called calibration 
load flow (CALF) based on FDEP has been proposed which is 
suitable for the real-time implementation of security fun- 
ctions such as int erconnection modelling and tie-line flow 
matching. The possibility of using FDLR or NETR in the 
place of FDEP with computational advantages should be examined. 

It should be possible to extend the two level technique 
proposed in Chapter ¥ for contingency evaluation. This work 
needs further investigation. 
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In the economic dispatch solution procedure using NETR 
proposed in Chapter IV, the line security constraints have not 
been considered. The inclusion of these additional constraints 
may pose difficulty. This needs to he examined further. 
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APPENDIX A 


ELEMENTS OP THE JACOBIAN MATRIX [J ] 


A. 1 POLAR COORDINATES 


The real and reactive power injections at the bus i in 

terms of the voltage angle and magnitude and the bus admittance 

matrix elements can be expressed as 
n 

p i " V i J ((G ij 005 9 ij + B ij sln u - 1) 

3-1 


n 


S i = V i 1 ((& ij sln 9 in - B ij 003 9 U )T 3 ) 


3=1 


ID ID 


ID D 


(A. 2) 


Then, the elements of the Jacobian [ J^J defined in equation 
(2.6) are as given below : 
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Vi fi _ 


3V i ~ G Ii V i + P i 


(A.6) 
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A. 2 BECTA1TGULAB COORDINATES 

The real and reactive power infections at the bus i in 
terms of the real and imaginary components of the voltage and 
"bus admittance matrix elements can he expressed as 
n 

I*i — I ^®i P f^ G i,j **" ^i ~ f-j) 


3=1 
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(A. 9) 


Q. = 7 (f. e. - e. f . )G. • - (e. e. + f. f.) B. . 
l * v l 3 i 2 13 i 3 1 3 13 

3=1 


(A. 10) 


The complex current flowing into the bus i is given by 
n 


* a i + 3 *i = (& lj + 3 B ij> (e j + 3 V 


( A. 1 1 ) 


Then, the elements of the Jacobian [J 1 defined in equation 
(2.11) are as given below : 
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APPENDIX B 


Mil B.IX MODIFICATION - TECHNIQUES 


B.1 INVERSION OF A MODIFIED MATRIX 


Let the changes in the original matrix [ B Q ] be 


represented as 


[A B] = 


i 1 b ii 


i f h . . 
J i D 1 


b. . ? 

ID i 


h . - ! 

DD ! 


(B.1) 


If = b^.. = b^ = b;._j = b, the new matrix fB-j] can be 
expressed as 


[B 1 J = [3 ] + b [ M] 5 { M] 


(B.2) 


where [M ] is a highly sparse row vector given by 
i D 

C 0 , 0,0 ... 1 , 0 , 0,0 ... 1 , 0 , 0 , 0 ] 

-1 rn -1 

Then [B-j] = { f B^ + b[Mj x [M] } (B.3) 

Using Sherman -Morrison formula [13] and expanding the equation 
(B. 3), we get 

“1 “1 - ^ T -1 T . 

CB^ =CB 0 ) -bCUbCM]^ [ M j ) ([B 0 ] [M] [M|[B 0 ] ) 

= - (£ +[ M][B q ] 1 [M?) 1 ([B o ] 1 [M ftMjpy 1 ) 

(B.4) 



168 


let [X ] = £ B Q 1 IM] T (B.5) 

a = (~ + tMHX] ) _1 (B.6) 

Using (B.5) and (B.6) in equation (B.4), we get 
-1 -1 -1 

[B-, J = (B o ] - a CXI [M][B o ] (3.7) 

Let the modified matrix equation be expressed as 

[B-jHX-,] = £Z 1 } (B.8) 

-1 

i.e. the solution vector [ X-j ] = [ B^ ] [ 2-j } (B.9) 

-1 


Substituting the expression for [ B-j ] 

from (B.7) in (B.9), we 

obt ain 




C X-j] ={ X q 3 + [ AXl 

(B.IO) 


-1 


■where 

{ X o 1 = [ B o ] C Z„ 1 

(B.11) 

and 

[AX ] = - a[X] [ M] [X 0 J 

(B.12) 


Solution Procedure : 

The solution vector £ X-j 3 can be calculated using the 
factors of the original matrix in the following manner. 

Step 1 : Compute the vector [X ] and the scalar 'a’ from 
equations (B.5) and (B.6) 

Step 2 : Compute [ X Q ] by a repeat solution using the factors 
of [B o ] in equation (B.11) 
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Step 3 : Compute the correction vector [A X] from eQn. (B. 12) 
and hence the solution vector t X-j } f rom eqn. (B.lO). 

B. 2 CALCULATION OP SENSITIVITY FACTOR 

In Chapter II, the consideration of Q-limit violations 
on the voltage controlled buses using sensitivity factors in 
the method of FDLR and NKM load flow techniques has been 
discussed. Here, the calculation of these sensitivity 
factors following reference £ 29 3 is given. 


Each P-V bus is checked for its Q-limits. If there is 
a violation, then the concerned bus is switched to a P-Q bus. 
That means, the matrix tB"3 , in equation (2.29) is modified 


by an extra row and column as 


EB"1 


new 




Let the order of the original matrix be (n-1 ) x (n-1 


(B.13) 

If 


the new P-Q bus is represented as n, then the order of 
[B"J ne w be hxn. The factors of the symmetric matrix { B” ] 

are expressed as 


fB"} = [ur 


c u 3 


(B. 14 ) 


where CU 3 is the upper triangular matrix with unit diagonal 
elements and fD] is the diagonal matrix. The factors of 
[B !, 3 new can Be obtained by adding the nth column to { U ] and a 
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diagonal element to [ in the existing factors. These 
additional factors can be obtained from the factors of the 
original matrix in the following manner. 

Solution Procedure : 

Step 1 : Obtain the intermediate vector [E-j] by forward 

T 

operation on [ C^] and [ U] . 

Step 2 : Obtain the additional factors U-j n » U^* **• ^n-1 n 
by multiplying [D J and [P^] and U nn =1. 

Step 3 : Obtain the nth element of [ D ] from 

D nn - (0 2 - V 4/V 1 (B - 15) 

d=i 

The element D Qn as obtained from (B.15) is known as the 
sensitivity factor S used for Q-limits. 

The above procedure can be followed with slight modi- 
fications for use in equation (2.46). Here, the matrix fJ] 
is not symmetric. Hence, both the lower and upper triangular 
factors of the original matrix need to be stored and operated 
to obtain the additional column as well as row factors for 
the calculation of sensitivity factors. 

B.3 MATRIX REDUCTION AND UPDATING TECHNIQUES 

In Chapter II, NEPR load flow solution using the network 
reduction technique has been proposed and in Chapter III, it 
has been extended for contingency evaluation using the matrix 
updating algorithm. 
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Here, we present the essential steps of the matrix 
reduction and updating techniques used in NETR load flow and 
contingency evaluation. The detailed derivation of these 
steps are given in reference [ 21 } 


B.3.1 Reduction process 

Consider the matrix of order nxn given hy 


[ I ] = 
1 n J 


m 


n-m 


; £ 


Y, 


(B.16 ) 


m 


n-m 


The matrix [ Y n 1 is reduced to l 1^] of order mxm hy applying 

T 

Gaussian elimination on [Y n J and eliminating ] and [Y£,J"by 
using an appropriate unit upper triangular matrix. 


Let 


CY 1 

n 


n-1 f Y 


an 


£ 

cn 


cn 


Y n „ i 
nn = 


(B.17) 


n-1 


and 


tv 


A i 

u n ! 


\ ' 


Y 

cn 

nn 


(B.18) 


1 
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Then 


C U Y U T ] 
n n n 


I 

where = Y^ - -“-T 00 

nn 


n-1 


i 

H ■ 


0 


t 

i 

I 

* 

I 

1 


0 


an 


(B.19) 


i.e. nth column and row of [Y ] are eliminated. 
Continuing the above procedure, we get 


[UY U T ] 
n 


I Y m 


0 


0 


l <~1 
1 a 

i 

I 

f 

\ 
» 


(B.20) 


where 


[ n ) = t u ffi+ i , u m+2 ... u n ] 


\ d m+1 


o ! 


[ D] = 


0 


n-1 


n i 


(B.21) 


> * 


U n 321(1 d m+1» d m+2» 


n : 


The values of U m+1 , U m+2 
obtained by Gaussian elimination procedure are stored for 
updating the matrix J for a given change in fY n j . 
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B.3.2 Matrix updating algorithm 


Let the change in matrix [Y Q ] be given as 


i * AY. . 

s IX 


AY. 


ID 


T 

- 


[A Y n ] = [^1 = j I 1 = [C 11 ] 


D i AY.. t AY. . : 

J 1 tl t) i 

l * -J 

i 3 (2x2) 


Let 


{ X n ) be a highly sparse matrix of order nx2 


(B, 22) 


5 0 0 1 

J a 

' 0 of 


! (n-1 )x2 


< • 

i 

j 1 

s 0 
: 0 
i 0 


0 

0 


ith row 


1 | jth row 


0 


(nx2) 


{ (1x2) 


(B. 23) 


and[h Q lbe a 2x1 vector to be used as a temporary storage. 
The main steps of the algorithm are as follows : 


1. 

2 . 

3. 

4. 


If i < 3 < n » then set k = 3; otherwise set k = n 

ih y = iz k j 


D k - d k t f^o 5 


lx 


■k-1 


J = C x£ 1 + t uj l x£ 1 
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5. CD k l =U t l-tI W l[^/l) k 

6. [O*- 1 ] = tS k ) - [h^][hy/D k 

k = k - 1 


If k > m, go to Step (2); otherwise, go to step (7). 

T 

7. Obtain the change IaTJ = [T* C m ] and add the 
change to [Yj . 



APPENDIX C 


COHVERGMCE THEORM 

Ihe convergence theorem used in Chapter VI is as follows : 

Let E "be continuously differentiable and JA(x)l be non- 
singular matrix with elements continuous in the ball 

01 = { x e IR N \ J [x-x 0 j| < % } 

I 

Let us define 

Contraction constant = 0 = max [| D <*> (£) || 

c £e$ 

Suppose 

1) 0 o < 1 

ii) ||* (x°) - I°|| < (1-0 o )e 

Then , 

A 

1) there exists a .unique solution of E(x) = Z in u< 

2) the sequence {x°,x\ ... > generated by <J> will 
converge to the fixed point x s of $ in a 

3) - x s j|<^ _ C G {jx k - x^ 1 j| 

c 

where IR = IT-dim easional vector space* 

it* ii = max {j^i> i.e. , <n-norm. on I# is used. 

| |A 1 1 = max { l | a. . |} 

i 3 13 

and D$(?) denotes the derivative of 4> at £ 
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DATA FOR TEST SYSTEMS 

In this appendix, the data for IEEE 14 bus, 30 bus 
and 57 bus test systems and for a heavily loaded 28 bus system 
[381 are given. Base MV A = 100. The data is compiled in the 
following manner : 

1) Bus data 

2) Line data 

3) Transformer data 

4) Shunt capacitor data 

5) Voltage controlled or P-V bus data 

6) Generator data for optimal load flow problem 

D.1 IEEE 14 BUS TEST SYSTEM 

D 1.1 Bus d at a 


Bus No. Generation Load 



m 

MVAR 

MW 

MVAR 

1 

0.0 

0.0 

0.0 

0.0 

2 

40.0 

0.0 

21.7 

12.7 

3 

0.0 

0.0 

94.2 

19.0 

4 

0.0 

0.0 

47.8 

-3.9 

5 

0.0 

0.0 

7.6 

1.6 

6 

0.0 

0.0 

11.2 

7.5 


continued 



Bus Fo 





7 

0.0 

0.0 

0.0 

0.0 

8 

0.0 

0.0 

0.0 

0.0 

9 

0.0 

0.0 

29.5 

16.6 

10 

0.0 

0.0 

9.0 

5.8 

11 

0.0 

0.0 

3.5 

1.8 

12 

0.0 

0.0 

6.1 

1.6 

13 

0.0 

0.0 

13.5 

5.8 

14 

* ■ M W 

0.0 

0.0 

14.9 

5.0 


D J 

1.2 Line data 


w m . w 

Line No. 

Buses 

Line 

impedance 

Half line 




charging 

Brom To 

R in p.u. 

1 in p.u. 

susceptance 





p.u. 

1 

1 2 

0.01938 

0.05917 

0.02640 

2 

2 3 

0.04699 

0.19797 

0.02190 

3 

2 4 

0.05811 

0.17632 

0.01870 

4 

1 5 

0.05403 

0.22304 

0.02460 

5 

2 5 

0.05695 

0.1 7388 

0.01700 

6 

3 4 

0.06701 

0.17103 

0.01730 

7 

4 5 

0.01335 

0.04211 

0.00640 


continued 
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Line Ho. 

Buses 

Brom To 

Line 

R in p.u. 

impedance 

X in p.u. 

Half line 

charging 

susceptance 

p.u. 

8 

5 

6 

0.0 

0.25202 

0.0 

9 

4 

7 

0.0 

0.20912 

0.0 

10 

7 

8 

0.0 

0.17615 

0.0 

11 

4 

9 

0.0 

0.55618 

0.0 

12 

7 

9 

0.0 

0.11001 

0.0 

13 

9 

10 

0.03181 

0.08450 

0.0 

14. 

6 

11 

0.09498 

0. 1 9890 

0.0 

15 

6 

12 

0.12191 

0.25581 

0.0 

16 

6 

13 

0.06615 

0.13027 

0.0 

17 

9 

14 

0.12711 

0.27038 

0.0 

18 

10 

11 

0.08205 

0.19207 

0.0 

19 

12 

13 

0.22092 

0.19988 

0.0 

20 

13 

14 

0.17093 

0.34802 

0.0 


D 1,3 Transformer data 


Transformer Buses 

Ho. From " To 


1 

2 


4 7 

4 9 


Tap 

setting 


3 


5 


6 


0.978 

0.969 

0.932 
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Bus No. 


D 1.4 Shunt capacitor data 


Susceptance p.u. 


9 


0.190 


D 1.5 Voltage controlled bus data (P-V buses) 


Bus No. Voltage 

magnitude p.u. 


Reactive^ power limits 
Minimum MVAR Maximum MVAR 


2 

1.045 

-40.0 

50.0 

3 

1 .010 

0.0 

40.0 

6 

1.070 

-6.0 

24.0 

8 

1.090 

-6.0 

24.0 


Slack bus voltage = 1.06/0 
D 1.6 Generator data 


Bus No. 

'min ' 

‘ ,v p max " 


cp ax 

*G 

Cost 

coefficients 


mi 

m 

MVAR 

MVAR 

a oi 

oi 

u - 
01 

1 

50.0 

200.0 

-20.0 

100.0 

105.0 

2.45 

0.005 

2 

20.0 

100.0 

-40.0 

50.0 

44.4 

3.51 

0.005 

3 



-40.0 

60.0 




6 

20.0 

100.0 

- 6.0 

- 45.0 

40.6 

3.89 

0.005 

8 



- 6.0 

45.0 





continued 
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Lower limits on voltage magnitude at all the buses including 
the slack bus 1 =0.90 p.u. 

Upper limits on voltage magnitude at all the buses excepting 
the slack bus = 1.1 p.u.; for slack bus = 1.06 p.u. 


Cost 

of 

generation F^ 

“ ^ a oi + b oi P Gi 

+ c oi 

?Qj_) $ per 

D. 2 

IEEE 30 BUS TEST 

SYSTEM 





D 2 

. 1 Bus data 



Bus 

No. 

Generation 

1 ^ 

Load 



MW 

MVAR 

MW 

' MVJLR 

1 


0.0 

0.0 

0.0 

0.0 

2 


40.0 

0.0 

21.7 

12.7 

3 


0.0 

0.0 

2.4 

1.2 

4 


0.0 

o 

• 

o 

7.6 

1.6 

5 


0.0 

0.0 

94.2 

19.0 

6 


0.0 

0.0 

0.0 

0.0 

7 


0.0 

0.0 

22.8 

10.9 

S 


0.0 

0.0 

30.0 

30.0 

9 


0.0 

0.0 

0.0 

0.0 

10 


0.0 

0.0 

5.8 

2.0 

11 


0.0 

0.0 

0.0 

0.0 

12 


0.0 

0.0 

11.2 

7.5 

13 


0.0 

0.0 

0.0 

0.0 


continued 
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Lower limits on voltage magnitude at all the buses including 
the slack bus 1 =0.90 p.u. 

Upper limits on voltage magnitude at all the buses excepting 
the slack bus = 1.1 p.u,; for slack bus = 1.06 p.u. 

2 

Cost of generation Ih = (a Qi + b Qi Pq^ + c Qi P &i ) per hr. 



_ V 

11 ' '• 



L.2 IEEE 

30 BUS 

TEST SYSTEM 





D 2.1 Bus data 



Bus Ho. 


Generation 

k r 

Load 


MW 

WAR 

MW 

WAR 

1 

0.0 

0.0 

0.0 

0.0 

2 

40.0 

0.0 

21.7 

12.7 

3 

0.0 

0.0 

2.4 

1.2 

4 

0.0 

0.0 

7.6 

1.6 

5 

0.0 

0.0 

94.2 

19.0 

6 

0.0 

0.0 

0.0 

0.0 

7 

0.0 

0.0 

22.8 

10.9 

8 

0.0 

0.0 

30.0 

30.0 

9 

0.0 

0.0 

0.0 

0.0 

10 

0.0 

0.0 

5.8 

2.0 

11 

0.0 

0.0 

0.0 

0.0 

12 

0.0 

0.0 

11.2 

7.5 

13 

0.0 

0.0 

0.0 

0.0 


cont inued 
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Bus Ho. Generation Load 


,, , 1 

MW 

mm 

MW 

mm 

14 

0.0 

0.0 

6.2 

1.6 

15 

0.0 

0.0 

8.2 

2.5 

16 

0.0 

0.0 

3.5 

1.8 

17 

0.0 

0.0 

9.0 

5.8 

18 

0.0 

0.0 

3.2 

0.9 

19 

0.0 

0.0 

9.5 

3.4 

20 

0.0 

0.0 

2.2 

0.7 

21 

0.0 

0.0 

17.5 

11.2 

22 

0.0 

0.0 

0.0 

0.0 

23 

0.0 

0.0 

3.2 

1.6 

24 

0.0 

0.0 

8.7 

6.7 

25 

0.0 

0.0 

0.0 

0.0 

26 

0.0 

0.0 

3.5 

2.3 

27 

o 

i 

o 

0.0 

0.0 

0.0 

28 

o 

• 

o 

0.0 

0.0 

0.0 

29 

0.0 

0.0 

2.4 

0.9 

30 

0.0 

0.0 

10.6 

1.9 
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D 2.2 line data 


Line No. 

Buses 
from' "To 

Line impedance 

R in p.u. v X in p.u. 

Half line 

charging 

susceptance 

p.u. 

1 

1 

2 

0.0192 

0.0575 

0.0132 

2 

1 

3 

0.0452 

0.1852 

0.0102 

3 

2 

4 

0.0570 

0.1737 

0.0092 

4 

3 

4 

0.0132 

0.0379 

0.0021 

5 

2 

5 

0.0472 

0.1083 

0.0104 

6 

2 

6 

0.0581 

0.1763 

0.0093 

7 

4 

6 

0.0119 

0.0414 

0.0023 

8 

5 

7 

0.0460 

0.1160 

0.0051 

9 

6 

7 

0.0267 

0.0820 

0.0042 

10 

6 

8 

0.0120 

0.0420 

0.0022 

11 

6 

9 

0.0 

0.2080 

0.0 

12 

6 

10 

0.0 

0.5560 

0.0 

13 

9 

11 

0.0 

0.2080 

0.0 

14 

9 

10 

0.0 

0.1100 

0.0 

15 

4 

12 

0.0 

0.2560 

0.0 

16 

12 

13 

0.0 

0.1400 

0.0 

17 

12 

14 

0.1231 

0.2559 

0.0 

18 

12 

15 

0.0662 

0.-1304 

0.0 

19 

12 

16 

0.0945 

0. 1 987 

0.0 

20 

14 

15 

0.2210 

0.1997 

0.0 

21 

16 

17 

0.0824 

0.1932 

0.0 


continued 
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• ’ * 

- - • ‘ " 

Line Ho. 

Buses 
Prom ’ fo 

line 
R in p.u. 

impedance 
' X in p.u. 

Half line 

charging 

susceptance 

p.u. 

• -■ 



■ — ■ - - 

- ■ - - — «- 

— ----- • 

22 

15 

18 

0.1070 

0.2185 

0.0 

23 

18 

19 

0.0639 

0.1292 

0.0 

24 

19 

20 

0.0340 

0.0680 

0.0 

25 

10 

20 

0.0936 

0.2090 

0.0 

26 

10 

17 

0.0324 

0.0845 

0.0 

27 

10 

21 

0.0348 

0.0749 

0.0 

28 

10 

22 

0.0727 

0.1499 

0.0 

29 

21 

22 

0.0116 

0.0236 

0.0 

30 

15 

23 

0.1000 

0.2020 

0.0 

31 

22 

24 

0.1150 

0.1790 

0.0 

32 

23 

24 

0.1320 

0. 2700 

0.0 

33 

24 

25 

0.1885 

0.3292 

0.0 

34 

25 

26 

0.2544 

0.3800 

0.0 

35 

25 

27 

0.1093 

0. 2087 

0.0 

36 

28 

27 

0.0 

0.3960 

0.0 

37 

27 

29 

0.2198 

0.4153 

0.0 

38 

27 

30 

0.3202 

0.6027 

0.0 

39 

29 

30 

0.2399 

0.4533 

0.0 

40 

8 

28 

0.0636 

0.2000 

0.0107 

41 

6 

28 

0.0169 

0.0599 

0.0032 
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D 2.3 Transformer data 


Transformer Bus? 3 Tap setting 

Wo. Prom” ' ' To 


1 

6 

9 

0.978 

2 

6 

10 

0.969 

3 

4 

12 

0.932 

4 

28 

27 

0.968 


D 2.4 

Shunt 

capacitor data 

Bus Wo. 



Susceptance p 


10 0.19 

24 0.04 


I) 2.5 Voltage controlled bus data 


Bus Wo. 

Voltage magnitude 

p.u. 

2 

1.045 

5 

1.010 

8 

1.010 

11 

1.050 

13 

1.050 


Reactive power limits 
Minimum MVJLR Maximum MVJLR 


-20.0 

60.0 

-15.0 

62.5 

-15.0 

50.0 

-10.0 

40.0 

-15.0 

45.0 


Slack bus voltage = 1.05/0 
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D 2.6 Generator data 


Bus 

No. 

p min 

^G 

^max 
' ^G 

n min 

Q G' ' 

<r 

Cost 

coefficients 

m 

m 

M7AR 

mm 

a . 

01 

V 

c - 
01 

• 

- — w 

■•■■■■■ 


■ — • * • 

- " 

- - 

- • * ■ 

1 

50.0 

200.0 

-20.0 

250.0 

0.0 

2.00 

0.00375 

2 

20.0 

80.0 

-20.0 

100.0 

0.0 

1.75 

0.01750 

5 

15.0 

50.0 

-15.0 

80.0 

0.0 

1.00 

0.06250 

8 

10.0 

35.0 

-15.0 

60.0 

0,0 

3.25 

0.00834 

11 

10.0 

30.0 

-10.0 

50.0 

0.0 

3.00 

0.0250 

13 

12.0 

40.0 

-15.0 

60.0 

0.0 

3.00 

0.0250 

/ * A - 



t/ r atr 


... 

1 

■ r a 


Lower limits on voltage magnitude at all the buses including 
the slack bus 1 = 0.95 p.u. 

Upper limits on voltage magnitude at all the generator buses 
excluding the slack bus =1.1 p.u. 

Upper limits on voltage magnitude at all the load buses and 
slack bus = 1.05 p.u. 

2 

Cost of generation 3?. = (a . + b . P n . + c . P n . ) £ per hr. 

X 01 01 brl 01 brl 


D.3 IEEE 57 BUS TEST SYSTEM 

3) 3.1 Bus data 


Bus No. 


Generation 


Load 



MW 

- - m* ‘ ,r jm-- _ 4* m 

mm 

MW 

M7AR 

1 

0.0 

0,0 

55.0 

17.0 

2 

0.0 

0.0 

3.0 

88.0 


continued 
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Bus No. Generation load 



MW 

MVAR 

m 

MVAR 

3 

40.0 

0.0 

41.0 

21.0 

4 

0.0 

0.0 

0.0 

0.0 

5 

0.0 

o 

* 

o 

13.0 

4.0 

6 

0.0 

0.0 

75.0 

2.0 

7 

0.0 

0.0 

0.0 

0.0 

8 

450.0 

0.0 

150.0 

22.0 

9 

0.0 

0.0 

121.0 

26.0 

10 

0.0 

0.0 

5.0 

2.0 

11 

0.0 

0.0 

0.0 

0.0 

12 

310.0 

0.0 

377.0 

24.0 

13 

0.0 

0.0 

18.0 

2.3 

14 

0.0 

0.0 

10.5 

5.3 

15 

0.0 

0.0 

22.0 

5.0 

16 

0.0 

o 

• 

o 

43.0 

3.0 

17 

0.0 

0.0 

42.0 

8.0 

18 

0.0 

0.0 

27.2 

9.S 

19 

0.0 

0.0 

3.3 

0.6 

20 

0.0 

0.0 

2.3 

1.0 

21 

0.0 

0.0 

0.0 

0.0 

22 

0.0 

0.0 

0.0 

0.0 

23 

0.0 

0.0 

6.3 

2.1 

24 

o 

* 

o 

0.0 

0.0 

0.0 


continued 
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Bus Wo. 

Generation 

Load 



MW 

1 - ‘WAR 

MW 

WAR 

25 

0.0 

0.0 

6.3 

3.2 

26 

0.0 

0.0 

0.0 

0.0 

27 

0.0 

0.0 

9.3 

0.5 

28 

0.0 

0.0 

4.6 

2.3 

29 

0.0 

0.0 

17.0 

2.6 

30 

0.0 

0.0 

3.6 

1.8 

31 

0.0 

0.0 

5 .S 

2.9 

32 

0.0 

0.0 

1.6 

0.8 

33 

0.0 

0.0 

3.8 

1.9 

34 

0.0 

0.0 

o 

• 

o 

0.0 

35 

0.0 

0.0 

6.0 

3.0 

36 

0.0 

0.0 

0.0 

0.0 

37 

0.0 

0.0 

0.0 

0.0 

38 

0.0 

0.0 

14.0 

7.0 

39 

0.0 

0.0 

0.0 

0.0 

40 

0.0 

0.0 

0.0 

0.0 

41 

0.0 

0.0 

6.3 

3.0 

42 

0.0 

0.0 

7.1 

4.4 

43 

0.0 

0.0 

2.0 

1.0 

44 

0.0 

0.0 

12.0 

1.8 

45 

0.0 

0.0 

0.0 

0.0 

46 

0.0 

0.0 

0.0 

0.0 


continued 
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Bus Ho. 

m 

Generation 

mm 

MW 

load. 

MVAR 

47 

0.0 

0.0 

29.7 

11.6 

48 

0.0 

0.0 

0.0 

0.0 

49 

0.0 

0.0 

18.0 

8.5 

50 

0.0 

0.0 

21.0 

10.5 

51 

0.0 

0.0 

18.0 

5.3 

52 

0.0 

0.0 

4.9 

2.2 

53 

0.0 

0.0 

20.0 

10.0 

54 

0.0 

o 

• 

o 

4.1 

1.4 

55 

0.0 

0.0 

6.8 

3.4 

56 

0.0 

0.0 

7. S 

2.2 

57 

0.0 

0.0 

6.7 

2.0 


D 3.2 line data 


line Ho. 

Buses 
From To 

line 
R in p.u. 

impedance 

X in p.u. 

Half line 

charging 

susceptanee 

p.u. 

1 

1 

2 

0.0083 

0.0280 

0.0645 

2 

2 

3 

0.0298 

0.0850 

0.0409 

3 

3 

4 

0.0112 

0.0366 

0.0190 

4 

4 

5 

0.0625 

0.1320 

0.0129 

5 

4 

6 

0.0430 

0.1480 

0.0174 

6 

6 

7 

0.0200 

0.1020 

0.0138 


continued ... 
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line Wo. 

Buses 

Irom To 

Line 
R in p.u 

impedance 
. X in p'.u. 

Half line 

charging 

susceptance 

p.u. 

7 

6 

8 

0.0339 

0.1730 

0.0235 

8 

8 

9 

0.0099 

0.0505 

0.0274 

9 

9 

10 

0.0369 

0.1679 

0.0220 

10 

9 

11 

0.0258 

0.0848 

0.0109 

11 

9 

12 

0.0648 

0.2950 

0.0386 

12 

9 

13 

0.0481 

0.1580 

0.0203 

13 

13 

14 

0.0132 

0.0434 

0.0055 

14 

13 

15 

0.02©' 

0.0869 

0.0115 

13 

1 

15 

0.0178 

0.0910 

0.0494 

16 

1 

16 

0.0454 

0.206 

0.0273 

17 

1 

17 

0.0238 

0.108 

0.0143 

18 

3 

15 

0.0162 

0.053 

0.0272 

19 

4 

18 

0.0 

0.555 

0.0 

20 

4 

18 

0.0 

0.43 

0.0 

21 

5 

6 

0.0302 

0.0641 

0.0062 

22 

7 

8 

0.0139 

0.0712 

0.0097 

23 

10 

12 

0.0277 

0.1262 

0.0164 

24 

11 

13 

0.0223 

0.0732 

0.0094 

25 

12 

13 

0 . 017 S 

0.0580 

0.0302 

26 

12 

16 

0.0180 

0.0813 

0.0108 

27 

12 

■ 17 

0.0397 

0.1790 

0.0238 


continued ... 
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Line No. 

Prom 

Buses 

tDo 

Line 

R. in p.u. 

impedance 

X m p.u. 

Half line 
charging 
suscept anee 

p.u. 

28 

14 

15 

0.0171 

0.0547 

0.0074 

29 

is 

19 

0.461 

0.685 

0.0 

30 

19 

20 

0.283 

0.434 

0.0 

31 

20 

21 

0.0 

0.7767 

0.0 

32 

21 

22 

0.0736 

0.1170 

0.0 

33 

22 

23 

0.0099 

0.0152 

0.0 

34 

23 

24 

0.1660 

0.256 

0.0042 

35 

24 

25 

0.0 

1.tS2) 

0.0 

36 

24 

25 

0.0 

1 .£30 

0.0 

37 

24 

26 

0.0 

0.0473 

0.0 

38 

26 

27 

0.165 

0.254 

0.0 

39 

27 

28 

0.0618 

0.0954 

0.0 

40 

28 

29 

0.0418 

0.0587 

O 

« 

O 

41 

7 

29 

O 

• 

O 

0.0648 

0.0 

42 

25 

30 

0.135 

0.2020 

0.0 

43 

30 

31 

0.326 

0.4970 

0.0 

44 

31 

32 

0.507 

0.7550 

0.0 

45 

32 

33 

0.0392 

0.0360 

0.0 

46 

32 

34 

0.0 

0.9530 

0.0 

47 

34 

35 

0.052 

0.0780 

0.0016 

48 

35 

36 

0.043 

0. 0537 

0.0008 


continued .. 
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Line Ho. 

Prom 

Buses 

To 

Line 

R in”p.u 

impedance 
. X in p.u. 

Half line 

charging 

susceptance 

p.u. 

-- - 

■- . 


"■ -> > ■ ■ 


- - — 

49 

36 

37 

0.029 

0.0366 

0.0 

50 

37 

38 

0.0651 

0.1009 

0.0010 

51 

37 

39 

0.0239 

0.0379 

0.0 

52 

36 

40 

0.0300 

0.0466 

0.0 

53 

22 

38 

0.0192 

0.0295 

0.0 

54 

11 

41 

0.0 

0.7490 

0.0 

55 

41 

42 

0.2070 

0.3520 

0.0 

56 

41 

43 

0.0 

0.4120 

0.0 

57 

38 

44 

0.0289 

0.0585 

0.0010 

58 

15 

45 

o 

• 

o 

0.1042 

0.0 

59 

14 

46 

0.0 

0.0735 

0.0 

60 

46 

47 

0.0230 

0.0680 

0.0016 

61 

47 

48 

0.0182 

0.0233 

0.0 

62 

48 

49 

0.0834 

0. 1 290 

0.0024 

63 

49 

50 

0.0801 

0. 1 280 

0.0 

64 

50 

51 

0.1386 

0. 2200 

0 

• 

0 

65 

10 

51 

o 

• 

o 

0.0712 

0.0 

66 

13 

49 

0.0 

0.1910 

0.0 

67 

29 

52 

0.1442 

0.1870 

0.0 

68 

52 

53 

0.0762 

0.0984 

0.0 

69 

53 

54 

0.1878 

0.2320 

0.0 


continued 
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Line Ho. 

From 

Buses 

To 

Line 
R in plu. 

impedance 

X in p'.u. 

Half line 
charging 
susceptance 
p.u. 

70 

54 

55 

0, 1 732 

0.2265 

0.0 

71 

11 

43 

0.0 

0.1530 

o 

i 

o 

72 

44 

45 

0.0624 

0.1242 

0.0020 

73 

40 

56 

0.0 

1.1950 

0.0 

74 

56 

41 

0.553 

0.5490 

0.0 

75 

56 

42 

0.2125 

0.3540 

0.0 

76 

39 

57 

0.0 

1.3550 

0.0 

77 

57 

56 

0.1740 

0.2600 

0.0 

78 

38 

49 

0.1150 

0.1770 

0.003 

79 

38 

48 

0.0312 

0.0482 

0.0 

SO 

9 

55 

0.0 

0.1205 

0.0 


D 3.3 Transformer data 


Transformer 

Buses 


Tap setting 

Ho. 

From" ’ 

To 


1 

4 

18 

0.970 

2 

4 

18 

0.978 ' 

3 

7 

29 

0.967 

4 

9 

55 

0.940 

5 

10 

51 

0.930 


continued . . 



Transformer Buses Tap setting 


No. 

Prom 

To 


6 

11 

41 

0.955 

7 

11 

43 

0.958 

8 

13 

49 

0.895 

9 

14 

46 

0.900 

10 

15 

45 

0.955 

11 

21 

20 

1.043 

12 

24 

25 

1.000 

13 

24 

25 

1 .000 

14 

24 

26 

1.043 

15 

34 

32 

0.975 

16 

39 

57 

0.980 

17 

40 

56 

0.958 


Bus No. 

18 

25 

53 


D 5.4 Shunt capacitor data 

Susceptance p.u 

0.100 

0.059 

0.063 
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D 3.5 Voltage controlled bus data (P.V. buses) 


Bus 

No. Voltage 

magnitude 

p.u. 

Reactive jo ow er _ limit s 

Min imu nf M VAR Maximum'’ MV AR 

2 


1.010 


-17.0 


50.0 

3 


0.985 


-10.0 


60.0 

6 


0. 980 


- 8.0 


25.0 

8 


1.005 


-1 40.0 


200.0 

9 


0.980 


- 3.0 


9.0 

12 


1.015 


-50.0 


155.0 





Slack bus voltage = 1.04/0 




D 3.6 

Generator data 


Bus 

No. 

pmin 

MW 

:.i "• » • . . jr _ 

■pin ax 
r G 

m 

^In 

MVAR 

^nax 

mm. 

Cost 

a - 
01 

coefficients 

b . c . 

01 Ol 

1 

100.00 

500.00 

-50.00 

200.00 

0.0 

2.00 0.00375 

2 

100.00 

400.00 

-50.00 

100.00 

0.0 

1.75 0.0175 

3 

50.00 

300.00 

-50.00 

60.00 

0.0 

1.00 0.0225 

6 

50-00 

300.00 

-40.00 

65.00 

0.0 

1.00 0.0250 

8 

50.00 

300.00 

-30.00 

100.00 

0.0 

3.25 0.0085 

9 

50.00 

300.00 

-20.00 

50.00 

0.0 

2.00 0.0042 

12 

50.00 

300.00 

-50.00 

200.00 

0.0 

3.00 0.0250 


continued 
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Lower limits on voltage magnitude at 

all 

the buses 


including the slack bus 

= 0.90 p.u. 




Upper limits on voltage magnitude at 

all 

the buses 


excepting the slack bus 

= 1.1 p.u.; 

for slack bus = 

1.05 p.u 

Oost 

of generation lb = 

(a oi + *oi 1 

Gi + 

°oi P Gi> £ 

per hr. 

D.4 

28 BUS HEAVILY LOADED SYSTEM 





D 4. 

1 Bus data 




Bus 

Ho* Voltage 

Generation 

Load 


magnitude 

MW 

MV AH 

MW 

' MVAR 


p.u. 





1 

1.05 

0.0 

0.0 

0.0 

0.0 

2 

1.05 

152.0 

0.0 

0.0 

0.0 

3 

1.05 

390.0 

0.0 

0.0 

0.0 

4 

1.05 

49,0 

0.0 

0.0 

0.0 

5 

1.05 

440.0 

0.0 

0.0 

0.0 

6 

1.05 

191.0 

0.0 

0.0 

0.0 

7 

1.00 

0.0 

0.0 

145.0 

50.0 

8 

1.00 

0.0 

0.0 

265.0 

-22.0 

9 

1.00 

0.0 

0.0 

00 

• 

o 

26.0 

10 

1.00 

0.0 

0.0 

384.0 

-2.0 

11 

1.00 

0.0 

0.0 

374.0 

34.0 

12 

1.00 

0.0 

0.0 

129.0 

3.0 

13 

1.0 

0.0 

0.0 

135,0 

58.0 


continued 


• • 



Bus No. 

Voltage 

Generation 

load 


magnitude 

MW 

MVAE 

MW 

MVAR 


p.u. 





14 

1.00 

0.0 

0.0 

60.8 

7.0 

15 

1.00 

0.0 

0.0 

118.0 

-18.0 

16 

1 .00 

0,0 

0.0 

75.0 

44.0 

17 

1.00 

0.0 

0.0 

139.0 

52.0 

18 

1.00 

0.0 

0.0 

0.0 

0.0 

19 

1 .00 

0.0 

0.0 

0.0 

0.0 

20 

1.00 

0.0 

0.0 

0.0 

0.0 

21 

1.00 

0.0 

0.0 

0.0 

0.0 

22 

1 .00 

0.0 

0.0 

0.0 

0.0 

25 

* 

1.00 

0.0 

0.0 

0.0 

0.0 

24 

1.00 

0.0 

0.0 

0.0 

0.0 

25 

1.00 

0.0 

0.0 

0.0 

0.0 

26 

1.00 

0.0 

0.0 

0.0 

0.0 

27 

1.00 

0.0 

0.0 

0.0 

0.0 

28 

1.00 

0.0 

0.0 

0.0 

0.0 


D 4 

. 2 line 

data 



line 

Buses 

line, 

impedance 

Half 

line 

No. 

Prom"’ So 

R in p.u 

” X in p.u. 

charging 
suscept ance 





p.u. 


1 

7 28 

0.0140 

0.070 

0 . 

0 

2 

8 27 

0.0078 

0.039 

0 . 

0 


continued 



Line Bu ses Li ne i mped ance Half line 

Ho. From To R in p.u.'* X* in p.u. charging 

susceptance 

p.u. 


3 

9 

26 

0.0058 

0.029 

0.0 

4 

10 

25 

0.0054 

0.027 

0.0 

5 

11 

12 

0.0070 

0.035 

0.0 

6 

12 

22 

0.0062 

0.031 

0.0 

7 

12 

23 

0.0010 

0.005 

O 

• 

o 

8 

13 

23 

0.0196 

0.098 

0.0 

9 

14 

24 

0.0230 

0.140 

0.0 

10 

15 

21 

0.0118 

0.059 

0.0 

11 

16 

21 

0.02 94 

0.147 

0.0 

12 

17 

18 

0.0230 

0.140 

0.0 

13 

17 

4 

0.0-258 

0.129 

0.0 

14 

18 

20 

0.0220 

0.110 

0.0 

15 

18 

22 

0.0062 

0.031 

0.0 

16 

19 

20 

0.0060 

0.030 

0.0 

17 

19 

27 

0.0296 

0.143 

0.0 

18 

19 

23 

0.0242 

0.121 

0.0 

19 

19 

3 

0.0070 

0.035 

0.0 

20 

19 

1 

0.0200 

0.100 

0.0 

21 

20 

21 

0.0063 

0.034 

0.0 

22 

20 

26 

0.0033 

0.019 

0.0 

23 

20 

2 

0.0112 

0.056 

0.0 

continued 
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Line 

No. 

Buses 

Prom $0 

Line impedance 

R in p.u. ’ X in p.u. 

Half line 
charging 
suscept ance 

p.u. 

24 

21 

22 

0.0220 

0.110 

0.0 

25 

22 

23 

0.0054 

0.027 

0.0 

26 

22 

6 

0.0070 

0.035 

0.0 

27 

22 

1 

0.0054 

0.027 

0.0 

28 

23 

24 

0.0072 

0.036 

0.0 

29 

24 

5 

0.0104 

0.052 

0*0 

30 

25 

26 

0.0038 

0.019 

0.0 

31 

26 

27 

0.0076 

0.038 

0.0 

32 

26 

28 

0.0128 

0.064 

0.0 



D 

4.3 Shunt 

capacitor data 


Bus No 

• 



Susceptance p 

.u. 

7 




1.2069 


8 




1.0252 


9 




2.0804 


10 




1.6665 


11 




1.5422 


13 




1.0042 


14 




0.2080 


15 




0.4476 


16 




0.7768 


17 




0.5868 
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